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1.  INTRODUCl'ION 


It  is  nov  generally  accepted  that  the  behaviour  of  a  reactive  material 
needs  to  be  considered  in  terms  of  its  explosive,  its  thermal  and  its  mech¬ 
anical  properties.  Most  explosives  deform  before  initiation  and  the  details 
of  this  deformation  determine  how  energy  can  be  localised  to  give  "hot  spots". 
Once  heat  is  produced  locally  the  "hot  spot"  temperature  is  controlled  by 
the  heat  being  produced  and  the  flow  of  heat  away  from  the  "hot  spot".  To 
understand  this  process  it  is  itr,port.ant  to  have  the  reaction  parameters  of 
the  explosive  and  thermal  properties  such  as  specific  heat  and  conductivity. 

For  many  explosives,  particularly  new  cempositions ,  these  are  not  always 
known  accurately. 

The  research  on  this  contract  has  considered  the  mechanical  and  thermal 
properties  in  some  detail.  Studies  of  initiation  and  propagation  have  been 
greatly  aided  by  high-speed  photography.  Research  on  reaction  kinetics  and 
thermal  properties  has  used  mass  spectroscopy,  differential  scanning  calori¬ 
metry  (DSC),  thermogravimetric  analysis  (TG)  and  scanning  electron  microscopy. 
Mechanical  properties  experiments  include  measurements  of  hardness,  yield 
strength,  coefficient  of  friction  and  fracture  surface  energy. 

The  four  areas  of  work  described  in  this  report  have  either  been 
published  (sections  2,5,4)  or  are  ready  for  submission  (section  I). 

The  texts  of  these  papers  are  given  in  full. 

2 •  the  IGNITION  OF A  THIN  LAYER  OF  EXPLOSIVE  BY  IMPACT;  THE  EFFECT  OF  ADDED 
POLYMER  PARTICLES.  G.M.  Swallowe  and  J.E.  Field. 

2 . I  Introduction 

A  study  of  the  deformation  behaviour  of  a  thin  layer  of  material  when 
impacted  is  relevant  to  a  range  of  problems  including  the  sensitiveness  of 
explosives.  A  standard  test  procedure  for  assessing  the  hazard  involved  in 
handling  explosive  niaterials  is  to  imp.act  a  sample  (typically  20-50  mg)  with 
a  falling  weight.  The  sensitiveness  of  a  sample  is  usually  expressed  in 
terms  of  a  502  drop  height  to  cause  ignition  (i.e.  the  height  from  which  the 
weight  would  ignite  50%  ol  a  series  of  samples).  There  are  standard  methods 
for  obtaining  the  50%  height  (see,  for  example,  Dixon  and  Massey,  1957).  However 
although  the  test  has  been  used  for  many  decades  it  is  only  recently  that 
a  systematic  attempt  has  been  made  to  photograph  t.he  sample  during  its 
deformation  (Heavens,  1973,  Heavens  and  Field  1974).  A  combination  of 
high-speed  photography  combined  with  pressure  measuring  techniques  helped 
establish  the  physical  processes  occurring  during  ignition  and  propagation. 

In  general  it  is  thought  that  the  initiation  of  an  explosive  by  mech¬ 
anical  shock  is  thermal  in  origin,  although  some  workers  believe  that  a 
tribochemical  or  a  molecular  fracture  me  .-rh  in  ism  may  be  responsible  in  certain 
circumstances  (Taylor  and  Weale,  1932,  Ubbolohde,  1948,  Fox,  1970).  On  the 
basis  of  localized  thermal  energy  or  "lot  spots"  as  the  source  of  the  ex¬ 
plosion,  four  possible  mechanisms  have  been  envisaged  for  ignition  by  impact. 

(i)  Adiabatic  compression  of  crappec*  gas  spaces  (Bowden,  Mulcahy,  Vines 
and  Yoffe,  1947,  Chaudhri  and  Field,  1974 >. 

(ii)  Viscous  heating  of  material  rapidly  extruded  between  the  Impacting 
surfaces  (Kirich  and  Tabor  1947,  Bolkhoviwinov  and  Pokhil  1958)  or  by  capi- , 
llary  flow  between  grains  (Rideal  and  Robertson,  1948). 
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(iii)  Friction  between  the  impacting  surfaces  and/or  grit  particles 
and/dr  grains  of  the  material  (Bowden  and  Gurton,  1949). 

(iv)  Localised  adiabatic  deformation  of  the  thin  layer  upon  i^chanical 
failure (Af anas' ev  and  Bobolev, 1971 ,  Winter  and  Field,  1975). 

It  was  shown  quite  convincingly  (Bowden  et  al,  1947)  that  the  presence 
of  gas  bubbles  affected  the  impact  sensitiveness  of  nitroglycerine, and 
Bowden  maintained  repeatedly  (Bowden  and  Yoffe,  1949,  1952,  1958,  Bowden, 
1950,  1963)  that  mechanism  (i)  was  the  primary  cause  of  ignition,  not  only 
in  liquids  but  also  in  solid  explosives.  This  view  has  with  equal  persis¬ 
tence  been  rejected  by  tlic  Russian  school  of  thought  (Andreev,  Maurina  and 
Rusakova  1955,  Bolkhovitinov,  1959,  Atanas'ev  and  Bobolev, 1971 ) .  In  impact 
experiments  on  solid  explosive  materials  it  was  discovered  (Kholevo,  1946) 
that  the  sensitiveness  was  considembly  reduced  if  the  sample  was  prevented 
from  flowing.  Afanas'ev  and  Bobolev (196))  produced  evidence  from  strain- 
gauge  experiments  which  suggested  that  during  impact  the  sample  suddenly 
undergoes  a  type  of  fracture.  These  authors  proposed  a  mechanism  of  hot¬ 
spot  formation  by  the  release  of  energy  along  slip  surfaces  formed  at  the 
onset  of  mechanical  failure  (Afanas'ev  and  Bobolev* 1971 ,  Afanas'ev,  Bobolev, 
Kazarova  and  Karabanov  1972). 

The  high-speed  camera  work  of  Heavens  and  Field  (1974)  showed  that 
impacted  samples  in  the  drop-weight  test  may  undergo  bulk  plastic  flow, 
show  evidence  of  partial  fusion  and  even  (with  PETN)  melt  completely.  The 
flow  speed  during  these  processes  is  considerable  and  may  reach  a  few  100  m  s 
Ignition  occurred  at  a  small  number  of  local  hot-spots.  Strain  gauge  measure 
ments  showed  that  pressures  of  typically  0.5  -  1  GPa  (5-10  kbar)  were  ach¬ 
ieved.  If  a  sample  failed  by  plastic  flow  this  was  accompanied  by  a  sharp 
pressure  drop.  It  was  after  the  pressure  drop  that  the  high  velocity  flow 
phenomena  took  place  and  ignition  occurred.  Hot-spots  were  thought  to  be 
caused  either  by  local  obstructions  in  the  flow  or  by  the  rapid  closure  of 
gas  pockets  trapped  in  the  material. 


Recent  work  has  been  concerned  with  studying  the  effect  of  added  grit 
particles  in  greater  detail.  The  case  of  the  hard,  high  melting  point 
particle  is  reasonably  well  understood  (Bowden  and  Gurton  1949,  BoWden  and 
Yoffe  1952,  Ubbelholde  1948).  The  basic  idea  is  that  when  two  solids  rub  or 
impact  together  the  hot-spot  temperature  at  the  interface  is  determined  by 
the  solid  with  the  lowest  melting  point.  Properties  such  as  thermal  con¬ 
ductivity  and  hardness  are  important  but  only  as  second-order  effects.  In 
fact  by  chosing  grits  of  different  melting  point  and  measuring  the  impact 
sensitiveness  Bowden  and  Gurton  were  able  to  measure  hot-spot  ignition 
temperatures  for  a  range  of  explosives.  The  hot-spot  temperature  found 
by  these  workers  for  PETN  in  a  drop-weight  impact  situation  was  ca.  700  K 
for  hot-spots  of  micron  dimensions. 

It  is  now  known  that  certain  polymers  such  as  polycarbonate  (PC)  and 
polysulphone  (PS)  can  also  sensitise  explosives  in  the  drop-weight  test 
though  others  such  as  polypropylene  (PP)  do  .not  (Bean  1973).  At  first 
sight  this  result  appears  surprising  since  polyir.ers  arc  relatively  soft 
and  have  softening  temperatures  usually  well  below  the  700  K  mentioned 
above.  The  primary  object  of  our  new  work  was  to  try  to  understand  the 
mechanisms  by  which  some  polymers  give  increased  sensitiveness  and  why  others 
do  not.  At  the  start  it  appeared  ponsiblc  that  the  explanation  could  be  a 
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chemical  one  baaed  on  the  production  of  free  radical a  which  aided  the 
reaction  kinetics  of  the  explosive  or  a  mechanical  one  in  which  hot**8pots 
were  produced  by  the  deformation  of  the  polymer.  The  first  experiments 
%fere  concerned  with  possible  chemical  effects  and  used  techniques  such  as 
differential  scanning  calorimetry  (DSC)  and  themogravimetric  analysis  (TG) . 
The  mechanical  effects  were  studied  using  the  photographic  and  pressure 
measuring  techniques  described  earlier. 

2.2  Studies  of  Chemical  Effects 

Experimental 

These  experiments  involved  measuring  the  reaction  parameters  for  decosr* 
position  of  pure  explosive  and  explosive/polymer  mixtures  and  assessing 
any  differences.  The  apparatus  used  was  a  Stanton-Redcroft  TG**750  thermo** 
balance  which  allows  the  weight  of  the  sample  to  be  monitored  whilst  it  is 
heated  at  a  constant  rate,  thus  giving  a  weight-loss  versus  temperature 
curve.  The  samples  used  consisted  of'^'I  mg  lots  of  a  mixture  containing 
25X  (by  weight)  of  polymer  to  75X  explosive.  This  percentage  of  polyiser 
is  much  greater  than  likely  to  be  found  in  practice,  but  was  chosen  to 
amplify  any  interaction.  The  explosive  used  in  all  these  experiments  was 
pentaerythritol  tetranitrate  (PETN).  The  analysis  follows  the  method 
of  Hauser  and  Field  (1978  and  section  3  of  this  report)  and  so  is  only 
briefly  outlined  here. 

Decomposition  is  assumed  to  follow  an  equation  of  the  form 


—  -  -  A  exp(-E/Kr)w"  (I) 

with  w  the  fractional  residual  weight  at  time  t,  E  the  activation  energy, 

A  the  frequency  factor,  n  the  reaction  order  and  T  the  temperature.  Since 
the  heating  rate  H  is  constant,  one  may  write  dT/dt  H.  Substituting 
in  equation  1  and  taking  logs yields: 

■  It  *  “ 

Thus  a  plot  of  (In  dw/dT  -  n  In  w)  against  I/T  gives  a  straight  line  if 
n  has  been  chosen  correctly.  The  value  of  E  is  determined  from  the  slope 
of  the  line  and  A  from  its  intercept.  The  analysis  was  carried  out  on  the 
C8iid>ridge  IBM  370/165  computer.  The  programs  used  fit  a  cubic  spline  to 
the  origiaal  data  and  produce  differentials  which  are  then  tiaed  to  give  plots 
of  (li^  -  n  In  w)  against  l/T  for  values  of  n  from  0  to  2  in  steps  of  0.2. 
The  straightest  section  of  these  curves  for  each  stage  of  the  reaction 
^here  may  be  sxire  than  one)  was  then  selected  by  eye  and  the  appropriate 
values  of  n,  E  and  A  obtained.  The  method  is  illustrated  by  the  activation 
energy  plot  for  a  sample  of  pentaerythritol  tetranitrate  (PETM)  in  Figure  I . 
The  technique  was  assessed  in  two  ways.  Firstly  by  using  the  obtained 
reaction  paraiMters  and  simulating  a  decosiposition  curve  which  could  be 
competed  with  the  original.  Secondly,  the  dehydration  of  calcium  oxalate 
monohydrate  was  observed. 

CaC20^H20  CaC20^  ♦  H^O 


This  reaction  is  well  understood  and  the  reaction  paraneters  are  well  docu¬ 
mented  in  the  literature  (Gurrieri,  Siracusa,  Cali,  1974,  Freeman,  Carroll, 
1958).  The  results  obtained  with  the  technique  outlined  above  were  in  good 
agreement  with  the  literature  values  and  are  set  out  in  Table  I. 

Table  I ;  Dehydration  of  Calcium  Oxalate  Monobydrate 

E/kJ  mol  *  In  A/ln  e  ^ 

This  work  89.5  -  3  20.2 

Literature  88.6  21.1 


Results 

(a)  Thermal  decomposition  of  PETN/polymer  samples 

The  polymers  reported  to  sensitise  PETN  (Bean,  1973)  were  polycarbon¬ 
ate  (PC),  polysulphone  and>er  (PS)  and  to  a  lesser  extent  polypropylene  (PP). 
Results  obtained  in  this  work  confirmed  the  dramatic  effect  produced  by  PS 
and  PC  but  showed  very  little  increased  sensitivity  in  the  presence  of  FP  . 
Samples  of  these  polymers  supplied  by  AURE  were  ground  in  a  freezer  mill  to 
sizes  of  <76  ym  and  then  mixed  with  PETN  to  produce  experimental  samples. 
The  heating  rate  used  throughout  was  10  deg  min  '  and  the  purger  gas 
nitrogen  or  argon.  The  results  are  sunmarised  in  Table  2,  the  figures 
listed  being  averages  of  between  three  and  seven  different  experiments. 


Table  2;  Effect  of  Polymeric  Additives  on  the  Thermal 
Deconyosition  of  PETN 

Sample  E  /kJ  mol  *  E./kJ  mol”*  In  A  /In  s  *  In  A./ln  s  *  n  T/K 
8  d  s  _d _ _ 


PETN 

143 

5 

193  -  7 

31 

♦ 

4 

44 

i  6 

0.7 

478 

PETN+PC 

130 

+ 

10 

214  -30 

30 

+ 

6 

53 

ilo 

1.0 

476 

PETN+PP 

122 

15 

172  -15 

28 

+ 

4 

41 

-  6 

0.6 

477 

PETN+PS 

130 

+ 

10 

231  -20 

30 

4 

49 

1 

MS 

O 

1.0 

474 

In  each  case,  a  two  stage  reaction  was  observed  (Figure  1),  the  first 
stage  being  identified  with  the  sublimation  of  PETN  (C  ,A  )  and  the  second 
its  decomposition  (E.,A.).  The  values  obtained  for  the  aStivation  energies 
of  subliaution  Q43  KJ  mol“*)and  decomposition  (193  *  7  kJ  mol“l)  are  in 
good  agreement  with  the  literature  values,  Cundall,  Palmer  and  Wood, (1977) 
quote  150.4  kJ  mol'*)  for  the  sublimation  of  PETN  while  Roberston  (1948), 
Rogers  and  Horris  (1966)  Haycock  and  Vemeker(t970)  and  Ng,  Field  and  Hauser 
(1976)  give  values  for  the  decomposition  of  197,198,188  and  192  kJ  mol**^ 
respectively.  The  value  for  In  A^  of  43.5  is  also  in  good  agreement  with 
the  45.6  reported  by  Robertson  (1948). 

A  larger  spread  is  observed  in  the  values  of  in  the  PETN/additive 
samples  than  the  pure  PETN.  However,  although  some  activation  energy 
differences  are  observed,  there  are  no  significant  changes  or  consistent 
trends.  The  values  of  frequency  factor  are  much  less  accurate  Chan  those 
of  activation  energy  since  they  have  been  obtained  by  smasuring  the  intercept 


*  The  sensitising  effect  reported  by  Bean  for  PP  in  PETN  is  surprising 
and  we  query  the  effect.  It  is  contrary  to  our  results  and  in  fact  tean's 
results  for  PP  in  other  explosives  show  little  or  no  positive  effect.. 
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rather  than  the  slope.  Again , within  the  errors  of  measureotent, they  show  no 
significant  changes  from  the  values  for  pure  PETK.  If  the  additive  had 
significantly  affected  the  reaction  towards  greatly  increased  reaction  rate 
then  A  should  have  increased  and  E  decreased.  Table  2  shows  that  only  small 
changes  in  A  and  E  took  place ,  and  further  that  any  changes  of  A  and  E 
tended  to  cancel  (i.e.  as  A  increases,  E  also  increases  and  vice  versa). 

The  only  consistently  observable  feature  in  the  results  is  a  change 
in  the  value  of  n  from  0.7  to  l.O  (except  in  the  case  of  PP).  A  value  of 
n  ca.  0.7  implies  a  surface  reaction  (for  which  n  =  2/3)  and  n  «  1 .0 
is  the  value  one  would  expect  for  a  bulk  reaction.  The  explanation  for 
this  is  that  the  high  proportion  (25%)  of  filler  material  distributed 
throughout  the  sample  breaks  up  the  PETK  giving  it  a  greater  surface  area, 
thus  favouring  a  bulk  rather  than  a  purely  .surface  reaction.  The  result 
for  PP  ir  explained  by  the  fact  that  it  softens  and  melts  at  a  lower  temp¬ 
erature  than  the  other  polymers.  PETN/PP  samples  were  observed  to  give 
droplets  of  polymer  which  then  coalesced  to  form  one  drop.  This  effectively 
removed  the  dispersion  of  particles  which  gives  n  approaching  unity  leaving 
a  surface  reaction  with  n  ca.  2/3. 

(b)  Thermal  decomposition  of  PETN/Benzoyl  peroxide  samples 

The  results  of  the  proceeding  section  show  that  the  additives  tested 
have  very  little  effect  on  the  thermal  decomposition  of  PETN.  It  may, 
therefore  be  concluded  that  their  sensitising  effect  is  not  due  to  polymer 
decomposition  products  interacting  with  PETN.  However,  the  possibility 
still  exists  that,  when  polymers  fracture  or  flow  during  impact,  large  con¬ 
centrations  of  free  radicals  may  be  produced,  and  these  could  react  with 
PETN  and  cause  ignition.  Walker  and  Green  (1976)  have  proposed  that  the 
mechanical  production  of  free  radicals  is  important  in  initiation  and  have 
done  some  experiments  in  which  they  find  that  tetramethylamroonium  tribo- 
hydride  (a  free  radical  donor)  causes  increased  reaction  with  anmonium 
nitrate. 


In  order  to  test  this  hypothesis,  mixtures  of  benzoyl  peroxide  and 
PETN  were  used  in  a  series  of  TG  experiments.  Benzoyl  peroxide  is  an 
organic  oxide  which  thermally  breaks  down  in  the  temperature  range  310 
to  350  K  to  form  two  identical  free  radicals 
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It  is  widely  used  in  the  polymer  industry  as  a  ’starter’  material  for 
free  radical  polymerisation  as  it  provides  a  good  .source  of  easily  formed 
free  radicals. 


A  number  of  TG  runs  were  carried  out  on  both  benzoyl  peroxide  by 
itself  and  two  PETN/benzoyl  peroxide  mixtures.  The  mixtures  used  had  com¬ 
positions  of  75%  to  25%  and  93%  to  7%  by  weight  of  PETN/benzoyl  peroxide. 

Due  to  the  decomposition  of  the  benzoyl  peroxide,  the  only  activation 
energy  that  can  be  estimated  is  the  decongjosition  of  PETN,  the  sublimation 
stage  being  buried  in  the  benzoyl  peroxide  reaction.  Simulated  computet 
decompositions  of  these  mixtures  were  made  by  combining  75%  of  the  nianerical 
value  of  pure  PETN  run  with  25%  of  the  numcricc 1  value  of  a  pure  benzoyl 
peroxide  run  (and  similarly  for  the  93/7  mix}.  The  figures  thus  obtained 
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the  usual  mnn\ji.'r.  Comuarisod  of  tlie  li,A  and  n  values  ohtninod  in  real 
and  simulated  erperimont  s ,  Table  T.  reveals  that  viiat  appears  at  first  to 
bo  a  very  strouj;  chemical  i.Ueraction  is  duo  to  the  e.ombined  effects  of 
the  two  materials  acting  independently. 
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(c)  Effect  of  ullra-violct  lieht  oi>  the  tljenual  decomposition  of  PETN/grit 
samples 

As  3  further  tost  of  the.  free-radical  initiation  hypothesis,  PETN/ 
additive  samples  were  irradiattd  with  ultra-violet  light  as  they  v?ere  being 
heated  in  the  TG.  ll-V  is  frequently  efloct.ive  in  degrading  polymers  by 
breaking  chains  with  the  fciae.ition  of  free  radicals,  both  PP  and  PS  arc 
susceptible  to  phetodegradnt  ifxi  by  IJ-V  whereas  PC  is  degraded  relatively 
much  less.  It  would  therefore  be  expcerevl  that,  if  free  radicals  are 
important  in  sensitising  PK.IM,  both  FP  and  PS  in  combination  with  U-V  should 
show  an  appreciable  effect  and  PC  little  or  no  effect.  The  results  obtained 
.in  these  experiment s  are  set  out  in  Tabic  A. 
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Except  for  the  ratlier  anonialous  behaviour  •>!  PS  in  the  decvampf's  i  r  ion  stage 
(;,ee  later),  t  I'.e  results  iidicate  that  "ri)  '  small  iTfects  ate  produced. 


The  changes  in  the  sublimation  parameters  E  and  A  tend  to  cancel  in  all 
cases  and  for  the  decomposition  parameters  and  changes  again  effect¬ 
ively  compensate,  particularly  for  polycarbonate.  PP,  PS  and  pure  PETN 

all  show  decreases  in  E^j  (Ed2  case  of  PS)  combined  with  increases 

in  Ad  (Ad2  with  PS).  Thus  the  reaction  is  helped,  albeit  a  little,  by 
the  U-V  radiation.  This  is  illustrated  by  the  slight  but  consistent, 
reduction  in  the  reaction  peak  temperatures  Tp. 

PS  shows  the  most  significant  changes  in  that  an  extra  reaction  stage 
appears.  This  result  was  found  to  be  very  repeatable,  though  it  is  not 
fully  understood.  Although  Edj  (for  the  new  stage)  is  very  high,  the  value 
of  Ad  I  more  than  compensates  for  this,  and  Tp  is  reduced  by  12  K.  The 
initial  and  final  stages  of  the  reaction  are  virtually  unchanged  by  the 
presence  of  this  new  stage,  which  only  has  a  small  (ca.  9  K)  range  of  activity. 

A  further  series  of  experiments  on  the  PETN/PS  system  was  carried  out 
in  which  the  U-V  light  was  switched  on  or  off  during  the  course  of  a  run 
and  the  effect  on  the  reaction  parameters  noted.  These  experiments  show 
that  the  stage  of  reaction  at  which  the  U-V  is  present  is  quite  important. 

If  the  U-V  is  switched  on  at  any  stage  before  a  third  of  the  material  present 
has  decomposed,  then  a  three  stage  reaction  is  observed,  otherwise  no  effect 
is  observed.  When  U-V  has  been  present  for  the  initial  stages  of  the  reaction 
(even  if  only  one  eighth  of  the  explosive  has  decomposed)  extinguishing  the 
light  has  no  effect  and  the  reaction  proceeds  in  three  stages,  as  if  the 
U-V  were  present  throughout. 

These  results  indicate  that  the  U-V  is  effective  in  breaking  down  the 
PS  to  produce  free  radicals  which  then  influence  the  decomposition  of  PETN. 

If,  however,  PETN  decomposition  has  proceeded  to  a  significant  extent  before 
the  U-V  is  present,  the  irradiation  has  little  effect.  This  is  probably 
because  the  free  radical  concentration  produced  by  the  PETN  alone  is  already 
much  greater  than  that  which  the  lamp  caiproduce.  The  observation  that 
extinguishing  the  lamp  early  in  the  reaction  does  not  stop  the  three  stage 
process  indicates  that  the  major  influence  is  the  concentration  of  free 
radicals  in  the  initial  stages  of  decomposition.  The  influence  of  U-V  may 
thus  be  summarised  as  producing  such  a  concentration  of  radicals  in  the 
early  stages  of  PETN  decomposition  that  the  decomposition  process  is  speeded 
up  and  the  reaction  parameters  (and  presumably  the  details  of  the  mechanism) 
are  changed. 

That  PC  is  not  influenced  in  this  way  is  not  surprising  since,  in 
comparison  with  other  polymers,  polycarbonates  are  not  degraded  to  any 
large  extent  by  U-V  (Ranby  and  Rabak,  1975).  It  is,  however,  surprising 
that  if  such  an  effect  should  be  produced  with  PS,PP  remains  ineffecitve. 

Both  PS  and  PP  are  degraded  to  a  large  extent  by  U-V  and  both  produce 
large  concentrations  of  radicals  when  degraded.  It  suggests  that  the 
sulphone  groups  have  a  specific  effect  on  at  least  one  of  the  stages  of  the  PETN 
decomposition.  Support  for  this  comes  from  the  work  of  Reich  (1973)  who  tested 
the  compatibility  of  a  range  of  polymers  with  various  explosives.  PS  is  one 
of  the  polymers  he  lists  as  incompatible  with  PETN  since  he  detected  an 
increase  in  the  value  of  the  reaction  parameter  n  in  the  present  of  PS. 


From  the  series  of  experiments  described  in  this  section,  it  may  be 
concluded  that  there  is  no  major  chemical  reaction  between  PETN  and  any 
of  the  polymer  fillers.  The  only  exception  to  this  is  PS  and  then  only  in 
the  extreme  conditions  of  having  a  high  intensity  of  U-V  shone  on  the  sample. 
However,  the  effect,  even  with  PS,  is  still  only  marginal  and  the  conclusion 
has  to  be  that  chemical  effects,  although  present,  are  only  of  secondary 
importance  compared  with  the  mechanical  ones  to  be  discussed  in  the  following 
section . 


2 .  1  Photographic  Observations  of  the  Impact  Behaviour  of  Polymers 
and  Explosives 

(a)  High-speed  photography  of  impacts  on  polymer  samples 

The  impact  behaviour  of  both  polymers  and  explosives  was  photographed 
using  a  CA  rotating  mirror  framing  camera  at  framing  intervals  of 

ca.  6.5  ps.  Impacts  were  observed  using  a  transparent  toughened  glass 
anvil  system  similar  to  that  of  Heavens  and  Field  (1974)  and  is  shown 
schematically  in  Figure  2.  The  weight  of  5.5  kg  was  dropped  by  an  elec¬ 
tromagnet  from  a  height  of  1  m,  the  fall  being  guided  by  three  cylindrical 
rods . 


The  behaviour  of  small  discs  of  polymer  ca.  2  mm  diameter  and  0.8  mm 
height  under  impact  is  illustrated  in  Figure  3.  Figure  3  shows  selected 
frames  from  an  impact  on  a  disc  of  PP.  A  trace  of  radial  expansion  versus 
time  is  given  in  Figure  5.  With  this  material,  there  are  no  rapid  changes 
in  the  expansion  rate  and  bulk  deformation  occurs  throughout.  Figure  4  is 
a  higher  magnification  sequence  for  an  initially  2  mm  diameter  disc  of  PS. 
Its  radial  expansion  is  also  included  in  Figure  5.  This  material  deforms 
initially  in  a  bulk  manner  but  eventually  fails  catastrophically  (see 
frame  4(c)  and  Figure  5).  PC  gives  similar  catastrophic  failure  to  PS. 
Microscopic  examination  of  the  samples  after  impact  confirms  that  catas¬ 
trophic  failure  of  PC  and  PS  is  associated  with  rapid  cracking  and  shearing 
in  localised  bands.  Figure  6  is  a  micrograph  of  a  deformed  PS  disc.  There 
is  a  network  of  fracture  throughout  the  sample  and  sets  of  fine  parallel 
shear  bands  (one  set  is  labelled  s).  An  enlarged  view  of  region  s,  but 
rotated  so  that  the  bands  are  horizontal  in  the  figure  is  given  in  6b.  It 
is  argued  later  that  it  is  these  regions  of  localised  deformation  (crack¬ 
ing  and  shear  banding)  which  give  rise  to  "hot-spots". 

These  observations  lead  to  the  conclusion  that  the  polymers  used  may 
be  classified  in  two  broad  groups:  (i)  Polymers  such  as  PC  and  PS  which  fail 
by  the  production  of  many  fast  moving  cracks  and  local  shear  bands.  This 
is  accompanied  by  rapid  radial  expansion.  (ii)  Polymers  such  as  PP  vdiich 
deform  plastically  at  quite  a  high  rate  but  do  so  by  bulk  deformation  and 
without  cracking  and/or  shear  banding. 

(b)  High-speed  photography  of  impacts  on  explosives 

Heavens  and  Field  (1974)  photographed  the  impact  of  a  wide  range 
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of  granular  explosives.  Powdered  PETN  was  observed  to  be  compacted  into 
a  pellet  and  then  to  undergo  severe  plastic  deformation.  The  material 
became  translucent  and  finally  completely  transparent  with  very  rapdid  f low 
(300  m  s~))  which  Heavens  and  Field  attributed  to  melting.  Ignition  occurred 
in  this  final  stage.  RDX  partially  fused  when  impacted  between  glass  anvils 
and  there  was  indirect  evidence  of  complete  sample  melting  of  RDX  between 
steel  anvils.  Pressure-time  curves  obtained  by  Heavens  and  Field  showed 
that  a  sharp  pressure  drop  occurred  when  the  sample  failed  plastically. 
Ignition,  in  samples  without  added  grit  particles,  always  occurred  after 
the  pressure  drop  and  when  the  explosive  was  flowing  rapidly.  In  the 
present  work  a  large  number  of  experiments  were  carried  out  on  pellets  of 
PETN  and  the  hshaviour  was  very  similar  to  that  described  by  Heavens 
and  Field  for  granular  material.  Figure  7  illustrates  a  typical 
result  in  which  a  25  mm  pellet  of  PETN,  pressed  to  5k  bar  and  of  diameter 
5  mm  is  impacted  by  a  5.5  kg  mass  from  a  height  of  I  m.  The  same  type  of 
behaviour  as  described  for  powders,  can  clearly  be  observed.  In  frames 
a  to  f  the  material  is  deformed  plastically  and  becomes  translucent.  Trans¬ 
parency  and  rapid  flow  (230  ms')  occur  in  frame  g  with  ignition  at  sites  A 
and  B  almost  simultaneously. 

The  results  obtained  with  a  large  number  of  such  sequences  show  that 
ignition  only  occurs  if  rapid  flow,  as  observed  in  frame  g,  takes  place. 
Although  rapid  flow  is  a  necessary  prerequisite  for  ignition  the  converse 
is  not  necessarily  true  and  many  failures  occur  without  a  resulting  ignition 
This  suggests  that  there  has  to  be  a  particular  region  in  the  flowing  material 
where  a  "hot-spot"  develops.  These  regions  could  be  associated  with  gas 
bubbles,  foreign  particles  or  enhanced  flow.  If  discontinuities  of  sufficient 
size  do  not  exist,  ignition  does  not  take  place.  Note  that  fast  reaction 
always  starts  at  localised  regions  (see  Figure  7  and  the  sequences  in 
Heavens  and  Field). 

(c)  High-speed  photography  of  impacts  on  explosive/polymer  samples 

A  number  of  high-speed  photographic  sequences  at  framing  rates  of  up 
to  2  x  10^  frames  per  second  have  been  taken  of  impacts  on  explosive/grit 
samples.  These  sequences  were  taken  in  transmitted  light  using  the  experi¬ 
mental  set-up  described  above.  The  first  point  to  emerge  from  these 
experiments  was  that, unlike  impacts  on  pure  PETN  where  ignition  never 
occurs  before  transparency  and  rapid  flow  of  the  sample, have  occurred, 
when  a  grit  is  present, ignition  can  take  place  before  the  explosive  has 
been  sufficiently  compacted  to  fail  plastically,  form  a  transparent  layer 
etc.  This  observation  implies  that  the  pressures  needed  to  cause  explosion 
are  much  smaller  than  with  a  pure  material. 

If  a  sample  ignites  while  it  is  still  opaque  it  means  that  small 
polymer  particles  within  the  pellet  can  not  be  seen,  and  thus  association 
of  the  initiation  sites  and  grit  particles  is  difficult.  To  overcome  this 
problem,  two  different  approaches  were  adopted.  In  the  first  instance, 
pressed  pellets  of  pure  explosive  were  surrounded  by  polymer  discs  and  the 
impact  phenomena  of  the  whole  array  photographed.  The  second  type  of  experi¬ 
ment  involved  pressing  pellets  of  explosive,  into  which  had  been  put  single 
polymer  discs,  effectively  modelling  a  single  large  polymer  particle  in  a 
sample  of  pure  explosive.  Examples  of  results  obtained  from  these  experi¬ 
ments  are  illustrated  in  Figures  8-11. 

Figure  8  shows  a  0.65  mm  high  PETN  pellet  surrounded  by  five  I  nm 
high  PS  discs  .  During  the  impact  the  PETN  disc  (dark  central  area)  breaks 
up  and  is  forced  into  close  contact  with  the  polymer.  The  polymer  discs 
deform  plastically  and  then  undergo  catastrophic  failure  as  described 
earlier.  In  frame  8b  all  the  discs  are  exhibiting  cracking  and  shear  band- 
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ing  and  by  frame  8c  have  failed  catatrophically.  In  frame  8d  we  can 
see  the  first  ignition  site  (labelled  I)  with  reaction  spreading  along 
two  paths  into  the  PETN.  A  second  site  (labelled  2)  has  developed  in 
frame  8e .  By  the  final  frame,  fast  reaction  has  spread  throughout  most 
of  the  PETN. 

An  example  from  the  second  type  of  experiment  is  shown  in  Figure  9. 

In  this  case  a  single  polymer  disc  of  PS  was  pressed  with  the  PETN 
powder  to  form  a  composite  disc.  The  selected  frame  shows  ignition  at 
three  points  (I  to  3)  on  the  polymer/explosive  interface.  The  polymer 
is  the  lighter  central  area,  and  the  darker  annulus  is  the  PETN.  In 
the  lower  part  the  PETN  is  becoming  ligher  and  this  is  due  to  the  onset 
of  fusion.  At  its  outer  edge  the  PETN  is  beginning  to  produce  fast  jets. 
The  mottled  appearance  of  the  polymer  is  partly  due  to  catastrophic  failure 
and  partly  due  to  a  very  thin  layer  of  trapped  PETN.  This  is  why  reaction 
also  spreads  inwards  as  well  as  outwards.  Frame  9b  shows  a  later  state 
where  fast  reaction  is  more  developed.  ThCjexpansion  velocity  of  the  PS 
disc  just  prior  to  ignition  was  ca.  35  m  s  .  The  catastrophic  failure 
occurred  between  frames.  Experiments  with  PC  showed  that  it  behaved 
in  a  very  similar  manner,  with  cracking  and  rapid  radial  expansion  followed 
by  ignition  at  the  interface. 

The  situation  with  PP  is  quite  different  and,  although  much  radial 
expansion  at  quite  high  velocities  is  observed,  the  polymer  does  not  crack 
or  shear  locally  and  ignition  does  not  occur. 

The  results  obtained  with  HMX/polymer  systems  are  similar  to  those 
obtained  with  PETN.  Figure  10  illustrates  an  impact  on  a  20  mg  sample 
containing  a  PP  disc  and  Figure  11  a  20  mg  sample  containing  a  PC  disc. 

In  both  cases  the  polymer  is  the  light  central  area  and  the  explosive  the 
darker  area.  Figure  10  shows  that  yield  and  flow  of  both  the  PP  and  HMX 
occur  but  no  ignition  results.  In  Figure  11  ignition  occurs  at  the  polymer/ 
HMX  interface  (frames  d  onwards)  after  catastrophic  failure  of  the  PC. 

The  findings  of  the  photographic  work  described  in  this  section  may 
be  summarised  as:  (i)  ignition  occurs  after  or  during  the  catastrophic 
failure,  cracking  and  fast  flow  of  the  polymer;  (ii)  the  initiation  site 
is  very  closely  associated  with  the  position  of  the  polymer  and  (iii)  of 
the  polymers  tested,  only  those  which  undergo  catastrophic  failure  are 
effective  in  promoting  ignition. 

2.1  Impact  Behaviour  in  the  Drop-Weight  Test 

(a)  Impacts  of  polymer  samples 

Samples  of  the  polymeric  materials  were  impacted  in  a  drop-weight 
machine  at  strain  rates  similar  to  those  observed  in  the  photographic 
work  (e  ca.  10  s  *)so  that  differences  in  their  mechanical  behaviour 
could  be  observed.  The  drop-weight  machine  used  was  similar  to  that  des¬ 
cribed  by  Heavens  and  Field  (1974)  and  is  illustrated  schematically  in 
Figure  12.  The  sample  was  placed  between  rollers  R1  and  R2  which  had  their 
edges  ground  down  to  give  a  sharply  defined  contact  area  of  diameter  5  mm. 
The  samples  used  were  discs  of  polymer  typically  1.5  mm  high  and  7  vm  dia¬ 
meter.  The  contact  area  was  lubricated  with  colloidal  graphite.  Two 
semiconductor  strain  gauges  were  mounted  in  series  diametrically  opposite 
to  each  other  on  R3  in  order  that  the  resistance  change  induced  in  the 
pair  of  gauges  during  impact  would  be  due  to  compressive  stresses  alone. 
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The  system  was  calibrated  statically  using  an  Instron  machine. 

As  the  sample  diameter  was  greater  than  the  contact  area,  the  area  of 
contact  during  impact  remained  constant  throughout  the  whole  process.  The 
pressure-time  traces  were  obtained  on  paper-tape  utilising  a  Data  Lab  DL922 
transient  recorder  and  analysed  on  a  HP9825  calculator  using  a  programme 
based  on  the  method  of  analysis  of  Afanas'ev  and  Bobolev  (1971).  The 
programme  numerically  integrates  the  force-time  trace  produced  by  the 
.equipment  to  yield  a  velotiy-time  trace  for  the  weight 

-1 

V(t)  -  V^-  M  F(t)  dt  (3) 


with  Vq  the  velocity  of  the  weight  at  impact  and  M  its  mass.  A  second 
integration  gives  the  displacement  of  the  weight  as  a  function  of  time  z(t). 

z(t)  -  V(t)  dt  (4) 

z(t)  is  the  compression  of  the  whole  system  (sample  and  anvils)  and  a  small 
amount  z  (t)  due  to  the  compression  of  the  system,  which  can  be  determined 
by  performing  impact  with  no  sample  present,  is  subtracted  from  z(t)  to 
yield  the  compression  of  the  sample. 

Zg(t)  ■  z(t)  -  z^(t)  (5) 

The  sample  strain  is  then  given  by 


c(t) 


In 


^  H 


(6) 


with  H  the  initial  height  of  the  sample  and  c(t)  the  strain  at'  time  t. 
The  stress  at  time  t  is  given  by  division  of  the  force-time  curve  by  the 
area  of  contact  A. 


a(t)  »  F(t)/A 

Elimination  of  t  between  equations  5  and  6  yields  the  stress-strain  curve. 

Results  obtained  with  PC  and  PP  are  illustrated  in  Figure  13  (a)  and 
(b).  The  behaviour  of  PS  was  similar  to  that  of  PC.  The  main  difference 
between  the  two  materials  is  the  evidence  of  catastrophic  failure  in  PC 
(as  shown  by  the  sharp  drop  at  e  ■  1.8)  and  its  absence  in  PP.  Both  mater¬ 
ials  show  evidence  of  strain  softening  followed  by  orientation  hardening, 
but  the  strains  which  PP  is  capable  of  are  far  greater  than  those  achieved 
by  PC.  The  onset  of  orientation  hardening  also  takes  place  much  earlier  in 
PC  at  a  strain  of  ca.  0.5  as  opposed  to  the  value  of  2.0  observed  with  PP. 

At  strains  of  >2  the  thickness  of  the  PP  layer  becomes  comparable  to  the 
correction  factor  for  anvil  deflection  and  thus  the  calculated  strain 
become  increasingly  unreliable. 

The  shape  of  the  stress-strain  curve  may  therefore  be  inaccurate  at  high  strains 
These  results  confirm  the  findings  of  the  photographic  work  in  showing 
a  marked  difference  in  high  strain  rate  impact  behaviour  between  sensitising 
and  non-sensitising  materials  and  again  point  to  catastrophic  failure  as 
being  the  process  responsible  for  a  polymer's  sensitising  action. 
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(b)  Impact  behaviour  of  explosive/polymer  samples 


The  sensitivity  of  pure  exp.losives  PETN  and  HMX  and  explosive/polymer 
samples  were  compared  by  determining  their  SOT.  height  on  the  drop-weight 
machine  outlined  above.  The  50%  height  is  the  height  from  which  the  hanraer, 
if  allowed  to  fall  freely,  would  have  a  50%  probability  of  causing  initiation. 
The  method  of  determination  was  based  on  the  "up  and  down"  method  which  is 
fully  described  in  Dixon  and  Massey  (1957). 


A  large  number  of  drop-weight  tests  were  carried  out  on  explosive/grit 
samples  to  determine  (i)  which  polymers  gave  rise  to  increased  sensitiveness, 
(ii)if  the  size  of  the  grit  particles  was  important ,(iii)  the  effect  of 
variation  of  the  percentage  weight  of  grit,  and  (iv)  if  a  polymer  which  sensi¬ 
tised  one  secondary  was  effective  in  sensitising  another.  All  the  experiments 
were  carried  out  using  twenty  25  -  1  mg  samples  of  powder  which  had  been 
pressed  statically  to  10  kbar  to  give  pellets  of  5  nan  diameter  and  approxi¬ 
mately  0.8  mm  height.  In  these  experiments,  three  grades  of  grit  were  used 
which  have  been  called  coarse,  fine  and  very  fine.  In  the  'very  fine'  grade, 
all  the  particles  had  been  through  a  76  pm  mesh,  the  'fine'  grade  indicates 
that  the  grit  was  used  as  obtained  from  the  freezer-mill  with  a  maximum 
particle  mass  of  ca.  20  Pg.  The  coarse  material  was  as-received  grit  from 
AWRE  with  particles  of  mass  up  to  1  mg  and  up  to  1  mm  long.  These  reuslts 
are  set  out  in  Table  5.  Unless  otherwise  stated,  all  the  samples  contained 
20%  by  weight  of  grit 

Table  5;  50%  Heights  for  Explosives  and  Explosive/Binder  Samples 


Sample 


50%  ht  (cm) 


PETN  (pure) 


21.8  -  1.0 


PETN  +  PC  (very  fine) 

PETN  +  PC  (very  fine)  repeat 


10.5  -  2.0 
10.7  -  0.6 


PETN  +  PS  (very  fine) 

PETN  +  PS  (fine) 

PETN  +  PS  (coarse) 

PETN  +  PS  (1  mm  high,  1mm  diam.) 


11.1  -  0.3 
13.0  ±  1 .5 

10.4  ±  0.6 

10.5  ±  3.0 


PETN  +  PMMA  (coarse) 
PETN  +  PMMA  (fine) 


21.0  t  1.6 

19.2  ±  0.8 


PETN  +  PP  (very  fine) 

PETN  +  PP  (very  fine)  10%  grit 


>8.8  t 
9.2  ± 


HMX  (pure) 


HMX  +  PC  (fine) 


16,0  -  1.5 


HMX  +  PS  (very  fine) 


15.1  -  1.0 


HMX  +  PMMA  (coarse) 


22.6  -  1 .0 


HMX  +  PP  (fine) 


24.4  -  0.8 


AgN^  (pure) 


19.3  =  0.8 


6.8  *  1.3 


AgN3  +  PC  (fine) 


From  the  results  set  out,  it  can  be  seen  that  the  size  of  the  grits 
used  does  not  have  any  great  effect  on  the  results  produced.  It  might  be 
expected  that,  if  the  effectiveness  of  the  grit  depended  on  its  ability  to 
disturb  the  flow  of  the  explosive  as  it  was  extruded  from  between  the 
anvils,  the  larger  particles  by  virtue  of  their  size  would  prove  more  effec¬ 
tive,  the  very  fine  particles  being  so  small  as  to  be  carried  along  in  the 
flow  of  the  explosive.  The  results  obtained  show  that  this  is  not  a  major 
factor . 

The  results  also  show  that  the  grits  are  not  specific  in  their  action: 

I'C  and  hS  being  equally  effective  in  sensitising  HMX  and  ’’ETN,  while  PP 
and  PMMA  only  slightly  sensitise  both  of  these  seconder'  explosives.  As 
noted  before  (see  for  example.  Heavens,  1973,  and  Heavens  and  Field,  1974) 
primary  explosives  such  as  lead  azide  Pb(N-).^  and  silver  azide  AgN^  are  not 
particularly  sensitive  in  the  drop-weight  test  provided  the  anvils  are  well- 
aligned.  However,  as  the  table  shows,  the  polymers  do  have  a  significant 
effect  on  the  sensitiveness  of  AgN^  and  in  a  consistent  way.  Finally,  in 
one  set  of  experiments  the  percentage  of  grit  was  reduced  to  10%  but  this 
caused  no  significant  change. 


2 . "  Thermal  Properties 

(a)  Mechanical  effects  and  thermal  properties 

The  mechanical  failure  properties  of  polymer  additives  appear  to  be 
of  major  importance  in  determining  their  sensitising  effect.  Those  mater¬ 
ials  which  undergo  bull  plastic  deformation  proving  ineffective  as  sensitisers 
and  those  which  fail  a  tastrophically  and  with  intense  local  deformation 
effective.  Thermal  pre oerties  areof  majoi  significance  in  determining 
the  failure  modes  of  pc l}nners  particularly  when  high  strain  rates  are  in¬ 
volved.  This  is  becauss  they  soften  at  relatively  low  temperatures  and  have 
low  densities,  specific  heats  and  thermal  conductivities.  Additionally, 
materials  can  desensitise  explosives  by  effectively  acting  as  heat  sinks 
and  quenching  hot  spots  before  they  can  fully  develop.  The  work  of  Bowers, 
Romans  and  Zisman  (1973)  for  example,  indicates  that  sensitivity  decreases 
with  increasing  specific  heat  of  additive,  those  additives  which  are  most 
effective  in  this  respect  being  materials  such  os "Superla  wax"or"Carbowax" 
which  have  low  melting  points  and  thus  absorb  their  heat  of  fusion. 

In  cases  where  thermal  properties  were  not  known  they  were  measured 
using  a  Perkin-Elmer  DSC-2  calorimeter. 

The  results  obtained  show  that  both  PC  and  PS  have  specific  heats 
which  range  from  1.25  kJ  kg”*  K  *  at  room  temperature  to  2.35  kJ  kg”'  K”' 
at  600  K,  the  increase  ^eing  reasonably  linear.  PC  has  a  melting  point 

of  420  to  430  K  which  oily  introduces  a  very  slight  increase  in j the  effective 

specific  heat,  giving  a  latent  heat  of  fusion  of  ca.  7.2  kJ  kg  .  PS  shows 
melting  over  the  range  470-500  K  with  a  heat  of  fusion  of  ca.  6.6  kJ  kg  '. 

PP  has  a  greater  specific  heat,  varying  from  1.92  kJ  kg”'  K”'  at  room  temp¬ 

erature  to  2.09  kJ  kg”'  K“'  at  450  K.  It  melted  over  the  range  390  - 
460  K  with  a  latent  hea:  of  fusion  of  ca.  64.0  kJ  kg”'  K“' .  PMMA  had  a 
relatively  constant  spe:ific  heat  of  2.9  kJ  kg”'  K”'  over  the  range  470  - 
630  K  where  a  large  decomposition  endotherm  was  produced.  From  the  results 
set  out  above,  it  could  be  predicted  that,  on  the  basis  of  thermal  properties 
alone,  the  order  for  ho  :-spot  quenching  for  the  polymers  investigated  would 
be  PP,  PMMA,  PC,  PS. 
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interface  to  be  by  an  infva-red  dclector.  Friction  camu'L  be  regarded 

as  just  a  surface  shearing  effect  since,  if  one  ef  the  irateriaTs  is  Imrder 
than  the  other,  tlie  nsperitie::  on  the  harder  surface  can  plough  out  the  softer 
material  to  an  appia-ci  able  dt'pl.li  below  the  surface,  it  is  thus  possible  for 
considerably  and  varied  dcforr-ial  ions  to  talce  place  in  a  frictional  situation, 
and  estimates  of  the  surface  temperatures  achieved  are  of  interest  both  to 
the  present  investigation  and  the  subject  of  polymer  friction  and  wear  in 
general. 

(a)  Apparatus  and  technique 

The  apparatus  constructed  to  perform  those  experiments  was  basically 
a  modification  of  that  described  by  Bov.’den  and  Thomas  (1954)  and  is  r,ho\\'n 
schematically  in  Figure  14.  It  consisted  of  a  50  mm  diameter  rotating  sapphire 
disc,  (infra-red  tranr,),iitting  to  ca.  5.3  pin)  against  v/liich  a  polymer  pin  vans 
loaded  liy  means  of  a  balance  arm.  The  area  of  the  interface  under  observation 
was  defined  by  a  0,8  mm  sl  it  mounted  directly  above  the.  rotating  disc.  The 
di.se.  rotated  at  a  frequency  of  ca.  56  Uz  giving  an  interfacial  velocity  of 
ca.  7  m  s"*.  Above,  this  disc  was  a  I’Ml-tA  chopper  wi.tli  radial  regions  in  the 
sequence  j  bl  acked-out ,  slit,  F>;ia\,  slit,  blacked  out  etc.  The  disc  v.ms  ro¬ 
tated  so  that  the  radi.ation  from  tiie  interface  was  cliopj'td  at  a  frequency 
of  10  kHz.  This  allot.’cd  temperature  variations  of  duration  1  ms  or  greater 
to  be  followed  (since,  for  any  one  temperature  determination,  a  complete 
bl acked-out  to  blackcd-out  sequence  must  be  obtained,  i.c.  5  'chops’).  After 
chopping,  the  radiation  was  delected  by  a  Mullard  liPY  36  infra-red  deleclor 
which  is  sensitive  over  the  r.aiige  1  -  5.5  pm  with  a  sharp  cut-off  at  5.5  pm. 
The  detector  v.’as  biased  to  give  the  best  .signal -t o-noi f.e  ratio.  The  detector 
output  vms  taken  Co  a  Tektronix  7A22  amplifier  and  the  amplifier  cutout  was 
fed  to  a  Datalab  DL9.'’2  transient  recorder.  Traces  v.’orc  then  punched  onto 
paper-t:’.pc  for  analysi.s  on  the  lJev;lett-Packrird  HP  9825  desl;-top  calculator 
and  plotter. 


The  purp.T.sc  of  the  ITv'A  filters  was  to  provide  .''n  absorption  ivand 
(ca.  3.5  pm)  v.’hi.eh  v.'ould  aiisorb  a  differing  propertien  of  the  signal  depend¬ 
ing  on  the  temperature  of  tic  source.  Thus,  by  calculating  the  ratios  I/T.. 
(with  I  the  inten.sity  llirougli  the  slit  and  I„  through  the  filter’)  for  sources 
of  various  knovT.  te::,;  eraturcs,  a  ca.iibraLier.  curve  of  1/lj.  ag.'iiirst  Ic.-ipc  raturc 
could  be  dr.’iv  n  and  estii.a-itcs  of  ter.riorater a  of  olin  r  sources  mad.e  using  the 
curve.  In  order  for’  the  calibration  to  bo  valid,  it  was  assumi:d  that  the 
sotneus  (cal  ibi’,1.  t  i  on  and  ;'j>ecir.ien)  were  In  h.  viug  ’  as  hi  ack- hod  i  er  .  Problems 
assoc’iatcd  v.’ith  di’.fe.rent  .sub.slances  having  differing  black-body  e’ni  ssivi  (  ies 
.are  c^lii.'inated  by  tli'O  use  of  the  r.rtio  I/lj.< 


(b )  Calibration 


The  theory  of  calibration  follows  the  treatment  of  Land  (l‘J64)  and 
Parker  and  Marshall  (19^8).  From  Planck's  radiation  law  we  have 

^  -  C  X'^  /  (e’'^XT  -  i)  (7) 

with  E;^  the  energy  emitted  by  the  black-body  at  wavelength  x  'Tod  temperatui'c 
T,  C  and  B  constants  and  the  eroissivity  of  the  surface  at  wavelength  X. 
This  can  be  approximated  to  Wien's  formula. 

E  =AGt:>  t  exp(-B/X  T)  (8) 

P  P  P 

with  E  the  ratliant  energy  received  by  the  detector, A  a  constant  dependent 
on  the  area  of  the  source,  C.  c  constant  dependent  on  the  geometry  of  the 
system,  t  the  transmission  coefficient  of  the  sapphire  and  X  the  effective 
wavelength  of  response  of  the  detector.  The  relation  betweeR  the  response 
of  the  detector  H  and  the  energy  falling  on  it  is  of  the  form 

H  =  E*"  (9) 

with  m  a  constant. 


Thus  from  equations  8  and  9: 
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In  H  ro(ln  (A  G  t)  ln<K,  X  )  “  B/X  T) 

X„  P  P 


(10) 


When  the  radiation  passes  through  the^PMUA  filter  the  optical  properties 
of  the  system  are  cluinged  (1,F.,X  e.nd  Xpcbange)  but  the  factors  A,  G  and 
t  remain  the  same.  Thus  if  H'  *'  is  the  reduced  output  on  passing  through 
•the  filter  we  may  write 


In  H'  »  m(ln(AGt)  +  Infe! 


X 


X’“^)  -  B/X'T) 
P  P 


From  equations  10  and  II 
H 


In 


H' 


m 


/,„(  li _ .‘i'l  .  I  il  -  i  U 

t  t  i  ^  ejJ 


<M) 


(12) 


This  relation  is  independent  of  tlic  source  area  and  geometrical  fractors  and  is 
independent  of  Che  absolute  values  of  the  omissivities.  It  depends  only  on 
the  differences  in  optical  properties  of  X  and  X'  (and  this  is  dependent  only 
on  the  properties  of  the  system  with  and  w?thouC  Ehe  filter).  Thus  a  plot 
of  In^,  .against  J.  .should  give  a  straight  line  and,  in  fact,  this  is  what 
is  found.  T 


Initially  caJibtation  vras  carried  out  using  a  2.5  mm  thick  PMMA  chopper 
disc  and  a  0.27  mm  diameter  platinum  wire  as  an  artificial  hot-spot.  The 
temperature  of  the  wire  was  varied  by  changing  the  IK  current  flowing  through 
it.  The  voltage  arrcis  a  central  portion  of  the  win’  w.s  monitored  via  tv.'o 
25.4  pm  diameter  platin'iin  wires  which  were  .spot-wi  ldeo  to  the  larger  wire  43.0  mr, 
.apart.  These  wires  were  thin  i  uough  not  to  affect  ll.c  tempt  ratairc  of  the  larger 
wire  which  was  consi ilcnd  to  be  uniform  over  the  48.0  run-,  long  middle  portion. 
From  the  knowledge  of  the  resist  ajKO  of  this  cent  ral  I'-trcion,  its  tt  mperaturo 
was  determined  using  standard  t  ihlcs.  The  ccntr.il  portion  was 


vit'Vl'cl  by  tlu'  ik'U’ctor  liuiu!;h  ;i  slit  of  vniith  1  inrr.  runintccl  porpciuli  riilar  ( c. 
tho  K’nyth  of  tho  ''.'iro.  Alllunii'l)  this  .’avc  n  goj-jd  linenr  rc'l  at  i  onsliip  oV'.  i' 
the  raiv’.e  770  -i;’7()  C  (Fip.nie  13a),  tb^-  ratios  lower  Le'njieratures 

bccai'.io  to  hiph  to  bo  Tneasxiroil  accurately  ami  errors  associ.ited  vv’ith 

tile  determination  of  current  ami  voltap.e  across  the  wire  beca:".e  sippii  ficant 
so  a  different  .'ippia-ach  was  a(l('pted  for  the  low  temperature  calibration 
of  the  system. 

Tile  thickness tho  PfC-lA  chopper  v/as  reduced  to  1  mm  in  an  attempt  to 
improve  the  ratio  and  a  variable  temperature  soldering  iron,  to  the  t  i )i  of 
which  was  attached  a  tlif'rnu'coupl e ,  w'as  used  as  the  source.  The  cnlibratioii 
produced  is  illustrated  in  Fip.ure  15(h).  A  good  1  ineay  re]  atlonsliip  i.s  ayain 
obtained  until  tiie  temperature  drops  below  ca.  500  F  (-  =  2.1  x  lO"*^)  when 

the  scatter  in  the  results  begins  to  get  cuiitc  considerable.  This  is  beca.in-.e 
the  ratio  has  again  incrcafU'd  to  values  where  the  aceuracy  in  its  de term  I na t ion 
is  reduced  by  noise.  At  such  lc>w  temperatures,  tlu'  powxir  emitted  by  a  stnin'e 
which  is  in  the  range  detectable  by  tho  detector  is  very  low  and  considerat i ons 
of  si  cnal -to-noTse  ratio  is  to  increase  the  size  of  the-  source.  A  simple  m.llv.d  . 
ot  iriprovinp,  the  signal  lo  noise  r.it  io  is  t('  increase  the  si/.e  of  the  sourm  . 

Tliis  ills  been  trieil  but  there  are  problems  assoeiate'd  with  dc'in;',  Sii  which  wi.!!  be 
described  below.  Furtlie^  reduction  in  tiie  thickness  of  the  I’lTlA  filter 
would  improve  Ihc'  ratio  yj,  ,  howi-ver  1  ima  is  the  smallest  commercially 
nv.-iilable  tliickiiess  and  it  wi'uld  be  nccersary  to  use  other  materials  to 
ir.iprovc  the  system  in  this  manner.  In  fact,  0.2  mm  thick  acetate  sheet  was 
tried  lint  w;is  found  to  he  unsuitable  as  the  position  of  its  absorpt  ion  bands 
was  sucli  that  vei'y  little  cb.aug.e  i.n  tho  value  of  II'  was  produced  for  a  wide 
range  (.'^OO-eOO  K)  in  T. 

The  calibration  v.ns  checluu!  by  measuring  the  radiation  received  frfim  a 
thernuii'ii  I  ei'  bulb  jilac.ed  benenlh  a  sapphiie  disc.  Tiie  bulb  was  Iteati'd  liy  a 
soldering  iron  and,  with  tlie  Iron  removed,  a  tempe.ratvire  reading  was  taken, 
and  the  result  obtained  from  the  infra-red  system  and  the  thermometer  scale 
compared.  For  a  tliermomcter  reading  of  K,the  infia-red  system  gave  a 
value  of  480  K.lbr  453  K  the  :>ystem  gave  438  K  and  for  the  low  temperatures 
of  383  k  the  sy.stem  gave  418  K,  It  can  thus  be  seen,  and  it  is  clear  from 
the  calibration  cur\’e,  that  the  determination  of  temperatures  below  ca.  500  K 
leads  t('  an  increasing  inaccuracy  in  the  result.  In  the  results  quoted  latir 
tho  error  in  the  risult  is  estimated  by  d''av2ng  an  envelope  which  encloses 
all  tlu  calibration  points  for  values  of  loss  than  2.1  x  10  (dotted 
in  Fip.uri.:  15h).  The  value  of  489  K  for  an  'actual'  temperature  of  488  K  may 
be  ratlur  fortuitously  close  and  all  estimates  of  temperature  will  have  an 
erren'  iif  the  order  of  *  10  K. 

Tile  prelitiiin  .ry  analysis  c'f  the  rei  ults  v;as  carrii'd  ^  ut  by  the  111’  9825 
cnleulati'r.  A  proi’.ramme  to  perform  this  analysis  was  written  which  changed  tbi' 
zero  level  of  the  traces  such  that  it  intersected  the  pcak-to-pcak  levels 
approx  i  mat  c’Ty  h.iir-WMV  (see  1  igpire  1b).  The  heights  of  the  peaks  were  then 
deduced  by  finding  the  halt-w.iy  points  h(  tw-on  the  intercepts  of  the  trace  and 
the  axis,  taking  .as  the  nni,i(>rical  value  of  the  peak  luig.'nt  an  average  of  S',\on 
))oints  .'rtiund  (.•.o'  iniliiding)  this  half-way  value.  The  trace  \.’as  then  phWiied 
out  ti  ;v  tber  will,  a  listiji,;  I'f  the  pc-a',;  ii.-!;,b!s  an.i  the  remaiiider  of  tlu'  aiia- 
ly.sis  (t  inding  I  lie  i  .at  io  was  done  bv  li.nidl.  Figure  lb  shows  three  lypioal 

resnlta.  (i)  is  tor  ,a  soliering  iron  at  '’7','  K  fii)  for  a  sr.iall  diaTi'eter  I’F  p,  i  a 

anil  (ii,!)for  a  larg.e  diar'eter  !'F  pin.  i’he  i  r'proveaienl  in  s  i  gn  al  -  ti  -noi  .se  i  a  I  io 

v.iili  a  larger  di  arx'ter  pin  is  clear  tri'ia  a  cony  .'w  i‘on  of  (iil^anc!  (iii)  la  la 

of  waicb  give  approx  iri.ate  I  y  Ihi-  s.an'c  nvi  n  i  auri  value  ot  tb.e  ratio 


ii' 


(c)  Results 


In  the  cxporinwntal  investigation,  polymer  pins  of  two  different  dia¬ 
meters  were  used:  snull  pins  of  diainetei  na.  1.3  mm  and  larger  pins  of  dia¬ 
meter  ca.  6  nEi'i.  The  smaller  diameter  pins  were  used  in  the  initial  exper i uienti; 
but  for  the  later  measurements  largo  pins  were  used  since  the  smaller  pins 
gave  a  poor  si gna  1 -t <v-r. ignal  noise  ratio  and  were  very  easily  deformed  due 
to  their  large  length/diameter  ratio.  The  main  problem  with  the  use  of  the 
larger  pins  was  that  eonsiderabte  ho.atiiig  of  thesipphiie  disc  took  place 
and  surf  ice  temporal  ires  of  up  to  360  K  were  rtcordeo  on  placing  a  thermo¬ 
couple  in  contact  wito  tiiC  disc  im-iediatel  v  (1-7  s)  .■ft>  r  doing  a  run. 

Due  to  the  con.siderabK-  raerinal  iPass  of  f'ne  disc,  its  rate  of  cooling  was 
quite  slow  (!C  -  20  d-giwes  per  minute!  and  sw  the  figure  quoted  above  gives 
a  good  indication  of  the  bulk  teniperature  reaclud  by  the  sapphire,  in  the 
course  of  an  exper i iis  r  t  . 

An  ei.perimerit  in  which  tl-.e  signal  level  v.-s  recorded  v.-ith  the  large  pin 
in  contact  with  the  disc  lud  again  i mmedia ! .  ■  ] y  alter  its  removal  indicated  that 
the  energtv  due  to  the  bulk  heeiing  of  tbe  disc  accouuicd  for  approxim-Htely 
one  thii'doj'  t  he  tote'  .  Such  ,i  ■'a>-pe  cont  ribut i on  from  a  lower 

temperature  .-.our<c  will  y.ive  a  mixed  sigtiai  viiich  wil  1  tend  to  reduce  the 
ateasured  blaek-hady  t otiipvrnt urc  ti  a  value  r.car  the  rnean.  Thus  the  readings 
c'btcineci  with  the  larger  pins  will  nrobsbl;'  ho  underirstinntes  of  the  interface 
tempera! 'iro .  In  the  :  iMo  of  results  giv.  i,  b  •  low  (Tahir  7)  the  maximum  black- 
b<)dy  tempoi at uro  rt-co- dod  by  tiie  svslem  for  iny  givci  polymer,  irrespective  of 
the  typr  of  pin  u-sed  is  quoted. 

The  rr-siiUs  wi  !  i  ,  (rf  '-..'•T  .  he  averages  over  tb.o  total  portion  of  the 

interface  area  vithi.  (he  dctcci.'-r’ ;•  'icld  of  view  aad  It  is  quite  likely  that 
s.mio  regions  (tlwcr-e  being  hignly  'b'.i c.nriod)  wi.’l  b.ivi'  temperatures  above  this 
avcr.igi;  valoo  and  otlu  .'s  will  h.ive  lover  v.aluc.s.  by  far  t!;o  largest  amount 
of  cxpoi'ir.cnta!  tine  spent  •n  tne  pc;  1  ymir:;  FV  -aiid  PC  re  t’nal  witit  an 

increased  anw  u  tt  of  ob-  i  rv.it  ion  lim*-  it  may  b,.*  possible  to  record  higher 
maxin'om  vnltH  s  of  ict.inrafue  fur  ibt  i  lirer'  c'thc!  rn.rteriai  s .  iris,  however,  in¬ 
teresting  tti.'.l  the  ronij  •  eaturvr  ->btainc;J  iollov  the  orucr  obtained  for  the 
fiensi  tivenes.-  .jf  th<.  p<' 1  aner /exp  1  ■ ■  vo  laixiuees.  It  is  al  .sc-  clear,  particularly 
in  the  ca'...-  of  FC ,  fi..ir  t  eiiipevat.ci  r  ;•  weli  abevo  the  coramerc  i  al  i  y  quoted 
sofleni'.ig  teripcratci  0  af  -iin  K  and  tl**-  r.bservcd  ’melt  ing’  temperature  (in 
previotis  work)  of  4.3!>  K  ar."  ubf  a  i n.tible .  In  th-.:'  c.ise  of  PF  however  this  wan  nor 
found  to  bo  tlu'  case  .ind  rc  tomperst  utei;  ab<r,L  the  commencement  of  melting 
at  IPO  K  wore  obr. air.t-d. 

Tabic  7:  Tcmpei^i^n'rj;  Aehj eve_d  Dviring  Polymer 


Pol  ymer  Maximuni  In  ter  fact  a]  Temperature/K 


♦ 

pp 

360  -  30 

IWA 

420  -  30 

+ 

PS 

430  -  .’0 

4' 

f’C 

490  -  10 

20. 


k 

k  / 

'  fk  * 

■n 


ICXiicun  '  I  i  or;  ol  I'.lic  polvivr’r  pin;:  nfl.i-r  I’hc  f  r  i  c  t  i  >'iia  I  f-'>:p(.‘r  iiiuao  s  ri-vc,.] 

Iha!  i  !i  l.lu'  ca.sc  oj’  I'C  .-iiiil  I’S  .‘.fvi’iro  (icliuri.it  U)!i  of  (hr  poTyiiicj;  haJ  booii 
prodiu  I'li  i.'iili  (.'v  i  lii'iic'-a.  <>i'  laiicli  s  I  rot  clr!  ii;;  .-.nd  slicair  i  i!-  of  molecular  layer.s. 

'J'lio  ri'  pin,  luv  ovor,  prost'iit  a  raacli  moi:o  roiiiuled  apyicaraiice  which  indicated 
Chat  lueJting  of  Llu'  ;;urface  layers  liad  taken  pliicc. 

Fi:,;ure  17  shows  the  overall  appearance  of  the  pins  after  the  cxpcrrimer.l  s 
(u)  is  a  PC  pin  and  the  many  layers  w'hich  have  sheared  off  and  built  up  on 
tl'e  trail  ii\''  edye  are  clearly  visible,  (b)  is  a  PP  pin  and  sliows  Llie  roundid. , 
mel  ted- 1  cokin.”  appe-arance  of  the  surface,  (c)  is  a  PMtf'v  pin  and  shovjs  the 
rou”bi.  iud  .siirfarc;  obtained  with  no  build  up  on  the  trailing  edge,.  All  the 
pins  in  Figure  17  were  of  diaineter  6  mm. 

M  i  c  ro.scopic  e;-:.’i..i  nation  eonrirmed  that  very  severe  deformation  of  PC  ha:! 
taken  place  with  clear  evidence  of  considerable  strain  in  the  frictional 
direction  (.Figure  18(a)).  PM'IA  (I'igure  lS(h))  presented  a  rather  'potjdery' 
appcaic.jce.  with  iiiany  p;vrticles  i' i  material  on  the  surface.  The  wearing  avciv 
of  tlu'  IMI.'i  pins  proceeded  i))  a  morelrittK;  manner  witli  small  particles  chipyu-;-! 
out  of  thr  surface  ai.d  with  no  a.ccuTmil a t i c'a  t'f  sheared  or  melted  layers  on  LiiC' 
trailing  eilge  cif  th.e  piit.  The  surface  of  I’P  on  the;  other  hand  was  relatively 
featureless;  as  niijuit  bo  expected  from  a  mcO  t;cd  and  refroxen  layer. 

Although  there  is  dear  ovicienco  that  PC  had  deforined  by  the  shearing  off 
of  la\’ers,  there  was  no  evidence  of  catas;  ropliic  failure  luic’ing  taken  place; 
the  shearing  cc'uld  be-  followed  by  tlu;  nakt'd  eye  and  took  place  over  a  time¬ 
scale  ot  ca.  1  s  i.k’r  a  6  mm  pin.  Kunethd  ess ,  ter.peraturcs  well  above  the 
comnicr  i  cal  softc'niii;,  t^nipi-'rat  ure  ■.■.'(mc  reec'rded  with  PC,  and,  v.’lu.n  aliowancf 
is  made  foi  ilie  fae.r  tl)at  these  wi-iu  averages  over  the  whole  siu'i  a.c.e,  it  does 
Seem  po-.i'ihle  that,  under  rendition;:  of  c->.trc;re  (and  pri'hably  locally  rapiil) 
deforr.'at,  Lon ,  very  liigh  t  oiviperaturcs  could  ho  achieved  by  i)o]yn:ers  t.ueh  as 
These  t erper.i tures  would  he  well  above  the  limits  set  by  the  old  ideas  of  the 
temp(.;ral  are  rise,  achieved  in  friction  being  limited  by  t!ic  melting  point  cif 
the  lower  meriting  i)oint  material  wliich  was  based  on  the  frictional  properties 
of  metals.  Any  increa;;c  in  melting  point  of  the  polymers  due  to  the  force 
applied  by  the  counterweight  would  be  negligible  since  the  maximum  woiglit 
used  v.'ith  the  6  mm  pins  was  2  kg  and  with  the  1.5  mm  pins  was  500  g  giving  a 
iTin;:jmum  stress  of  3  mP.-i,  which  corresponds  to  ca.  30  bar.  This  would  have  an 
insignificant  effect  on  the  melting  point,  variations  of  ca.  200  bar  being 
needed  to  cause  an  ai'preciable  (ca.  10  K)  change  (Zoller  1978).  It  does 
tht;referi’  appear  that  high  temiicratures  can  In;  reached  during  the  deformation 
of  poly.aers,  and,  in  the  design  of  any  system  containing  a  polvT.icr,  considerali 
must  be  givin  to  the  p-'S.s  ibi  1  i  ty  of  high  temperature  regions  during  defori;;  at  iou 

2.7  Coin  1  u;;  I.('ns  aprl  Iki  r-cu;'s  I  on 

A  c  or:h"  nut  i.in  o.'  hi  j’.h-spccd  j'hetograph  !  c  ,  pressure  measuring  and 
tlurr.,;!l  aiir.lNv,!;;  tcehnifuies  have  been  us^'d  to  study  tin,'  behaviour  and  inter¬ 
act  i.ri  ot  i;;pl  o.s  i  \  ai'd  ’p<jl  y,',  i  a  .  f.lu'^'.ical  in  t  iv  r  ,i  c  t  i  ('a  s  ,  for  the  imtcrial;. 
tesli.c’,  were  r.'iiiij  (o  be  ef  miiio;  :  ::'p.  r  t .  uc' ■  co’iipared  vi  t  li  mv'chanie.l  c  en.'- i  o'- 
eralo  u.  •  ir  ir-p;.  ;  si.;,:  i  t  i a  t  i  i  .  Pure  (Cplo.sivcs  v.'ei'e  found  to  ignite  onl  v 

after  r.'p  il  .iiu  u.  d'.’i’,  ion  of  the  .'..rpl-'  had  occurri'd.  Tn  the 

pn.ea  I:  ■  ,,r  ;  or,  ,  i  ;  i  i  to;  '.rit',  r  ,  ta.;,,  si  vt  I'e  ib  ■  t  e  t  i  on  was  not 

I"''  '  ■  '  ..:y  tor  i  n  :  1  i  or. .  bolb  ■  I'o  i  r  1  li  la  r.r.d  tin  rm'l  prc'pcrt'of  I'f  an 

p '  -  hi  '  '  i  'a  i  p  ;  I  in  d.  I  ,  ■ .  .ling  I  !  ■;  !  s  i  t  i 1  I  v  a  f  i  p.s '  <  .  A  iv  1 1  vra  r 

b  ■  '  p  '  b  1 1  a ;  1 . '  ■  p  ppI’  il  i  I ;  !i  '.-/it!:  I'v.-  sp,ai!:p'  ii.at  .md  ii:'-.'  ivat' 

'  !  t  u  o  ' ,  .'k  ,  di  ;  .  .  d  1  p  p  ■  i-  V-  ■ :  •  ■  r.p  1  .a,.  I  1  :  n'l  !  e  f  a  i  1  u  :  .  ,  v.'  i  11  '.u-  1  ;  \  p.'  i  y  to 


cause  sensitisation.  In  order  to  tost  this  hypothesis,  two  polymers  were 
chosen  on  the  basis  of  their  strength  and  thermal  capacity;  polyester  which 
has  high  strength  and  low  capacity  and  high  density  polytherfe  (HDPE)  v;hich 
has  a  low  strength  and  high  capacity.  They  were  mixed  with  PETN  and  tested 
in  the  drop-weight  machine  in  the  same  mao’.ier  as  the  samples  described  earl  icr. 
The  results  obtained,  together  v/ith  a  few  comparison  results  from  Table  5  are 
set  out  in  Table  8.  Thermal  properties  can  be  found  in  Table  6. 

Table  8:  50%  Weignts  for  Expl osive/Pclymor  Samples 


PETN  (pure 


PETN  +  HOPE 


PETN  +  PP 


PETN  +  Polyester 


PETN  +  PS 


50%  height  (cm) 


21.8  -  1.0 


24.8  -  1 .5 


18.8  -  0.8 


15.5  -  1.0 


n.l  -  0.3 


These  results  confirm  the  predictions  that  polyester  should  sensitise 

and  HDPE  should  not.  In  fact  it  appears  that  IIDPE  is  acting  as  a  desensitiscr . 

The  conclusion  reached,  therefore,  ir  that  polymers  which  are  most 
likely  to  .sensitise  are  those  sucti  as  PC  which  are  tough  below  their  glass 
transition  temperature  (a.s  opposed  to  PMMA.  which  behaves  in  a  much  more 
brittle  manner),  and  are  likely  to  fail  cat.istrophicully .  Those  least 
likely  to  sensitise  are  those  which  have  a  low  glass  transition  temperature 
(PP  or  HOPE)  and  which  deform  in  bulk  with  little  local  concentration  of 
energy.  Polyester  i.*^  somewhat  intermediate  in  behaviour  (and  sensiti,sing 
effect)  between  PC  and  PMtiA  as  it  it  tougher  than  PMMA  but  much  more  brittle 
than  PC,  and  so  fits  quite  well  into  the  ixpectcd  pattern. 


In  the  case  of  the 
important  property  (Kuwd 
tially  bccau.se  tliey  are 
relatively  l.irge  latent 
point.  In  contrast  th- 
a  range  of  temporat  i.rcs 
If  beat  is  proriured  by  loc 
However,  wlicn  softei-rng 
maintain  a  high  viscuf.it 
tribute  further  temper. it 
lo.s.ses  in  the  case  of  .i 
material.s  because  cf  (ho 
and  the  hot  spot  is  not 


former  the  melting  point  of  tiie  particle  is  the 
en  an.!  Curton  !h49).  That  titis  is  so  is  cssen- 
matcninls  with  well-defined  melting  points  and 
lieats;  the  hot  snot  is  generated  to  the  melting 
(lolymcric  m.iturial.s  which  sensitise,  soften  over 
and  h.'cve  low  latent  hears  and  thermal  conduc^tivitios. 
alisc'd  ratasJ  ruphic  failure  hot  spots  will  form, 
of  the  polymer  starts  the  material  will  still 
y.  Continuing  deformation  will  be  able  to  e.on- 
ure  rises  by  viscous  heating.  Thor  thermal 
polymer  are  very  low  comjiared  with  most  other 
low  latent  heat,  thermal  capaciry  and  conductivity, 
quenched  so  readily. 


The  friction  experiment  gave  evidence  of  temperalure.s  well  above  the 
.St)f  leu  i  iig  poiiii,  1 1;  this  experiment,  an  average  temperature  for  a  contact 
area  of  ca.  25  mr.:-  w.js  rwainlaincd.  I.ucal  Imt  spot  temperature  v/ould  be 
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expected  to  exceed  the  values  recorded.  In  other  experiments  in  the  laboratory, 
the  temperature  rise  at  tiie  tip  of  a  crack  piopajjat ing  in  a  polymer  wa.s  re¬ 
corded,  again  using  I  R  techniques  (Fuller,  Fox  and  Field  1973).  Two  impor¬ 
tant  results  emerged  from  this  work.  first  was  that  temperatures  well 

above  the  softening  points  were  observed.  Secondly, that  Che  fracture  surface 
energy  increased  with  crack  velocity.  This  latter  result  was 

explained  by  the  high  temperatures  causing  softening  of  the  material  at  the 
crack  tip  thus  requiring  more  energy  to  be  expended  in  crack  growth. 


It  is  clear  from  the  earlier  work  and  the  present  work  with  explosive/ 
polymer  samples  that  some  polymers  can  reach  very  high  local  temperatures 
during  their  deformation.  This  is  clearly  relevant  to  the  safe  handling 
of  explosives  since  polymeric  materials  may  be  used  for  packaging  or  as 
binding  materials  in  explosive  composites.  The  results  also  have  interest 
to  studies  of  the  fracture,  friction,  erosion  and  wear  of  polymers. 
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FIGUKK  CAVnmS 

1.  Acti.vaf.ion  en«'rf;y  plot  for  the  thcrm.il  docompositi on  of  PETK  for 
v,Trioiis  VvTIuos  of  react  ioo  ordfr  o.  'Ihe  n  vnlvit  vhich  give..s  the  best 
straiglit  lino  scr.rnunts  is  one  .qipropr  i  at  c  to  the  reaction  (here  0.8). 

Those  two  stops  (see  dashed  lino)  c<jrrospr)nding  to  sublimation,  stage  I, 
and  docomposi i ion  stage  11. 

2.  Experimental  arrangement  at  instant  of  impact,  W;  drop-weight;  G, 
glass  blocks;  M,  mirrors;  S,  sample.  The  tipper  glass  block  is 
attached  to  the  weight. 

3.  Impact  on  a  disc  of  PP.  Thi:  sample  deftirms  by  bulk  plastic  flow 
throughout.  Diamoter  of  field  ot  view  20  mm.  Frame  tirne/ps  ; 

a, ();  l).70;c,2fi6.  The  zero  tinges  given  in  this  and  subsequent  sequence.s 
do  mn  correspond  to  time  after  initial  contact. 

A.  Impact  of  a  disc  of  PS.  The  .s.imple  deform.s  initi.illy  by  bulk  plastic 
flow  (frames  a  na<'J  b)  hut  ther,  deforms  catastrophically  (frame  c) . 

Diameter  of  field  of  view,  2.1  mm.  Franu*  times/’is  :  a  ,0  ;b  ,2  I ;  c  ,  28 . 

5.  Diameter  versus  time  traces  for  impact. s  on  PP  and  PS.  Note  th.at 

the  diameter  of  the  PK  has  approximat-'ly  dt'ublcd  before  the  catastropic 
failure,  but  that  the  PP  di.se  ha.s  a  fin.a1  diameter  several  times  its 
initial  value. 

6.  Micrographs  of  a  doformed  PS  di.sr  .showing  lhat  the  catastrophic  failure 
is  a  result  of  fr.actnre  and  localised  she.ar.  '’'he  enlarged  view  in  (b) 
is  from  region  S  in  (.i).  The  view  has  been  rotated  so  that  the  shear 
bands  an.  horixont.al . 

7.  Impact  on  a  pellet  of  1'F.TN.  Franit*  a  shows  the  5  mta  diameter  pellet 
before  iraimcl.  P.apid  flow  .-nid  jetting  occur  as  for  powdered  layers. 

Fusion  of  tht  layer  starts  in  frame  d.  The  wliole  layer  is  transparent 
and  flowing  ra\>idly  by  f r.unc  g  when  the  first  to  two  ignition  sites 
forms.  The  circui.ii  feature  in  frame  h  i.s  an  artifact.  Diameter  of 
field  of  view  20  mm.  Frame  times/ps  :  .a ,0 ; b ,67 ;c , 22 1  ;d , 3C8 ;e , J89 ; 

f  ,A69;g,476;h,.'t82;i4?.9. 

8.  PF.TN'  pellet  surrounded  by  b  FS  discs  Sec  text  for  details.  Two  ig~ 
nition  sites  at  the  polymer/explosive  interface  are  labelled  in  frames 
d  and  e.  Diameter  of  field  of  view  20  min.  Irame  times/ys  :a,0; 

b ,  9 1 ; c 1 40 ; t , 16  1 ; c , 1 68 ; f , i 7b . 

9.  25  mg  disc  of  PETN  ronlaining  a  PS  disc.  Three  ignition  sites  are 

labelled.  Field  of  vuew  20  mro.  Frame  tintes/ys:  a,0;b,!4. 

10.  A  PP  disc  (light)  pressed  into  an  HMX  pellet.  The  first  frame  is  before 
impact.  Considerably  bulk  deformation  of  t'ne  PP  occurs  but  there  is 

no  ignition  of  the  explosive.  Field  cf  view  20  mm.  Frame  times/\is: 
a,0;b,9l ;c, 161 ;d , 2 1 4 ;e ,247 ; f ,402 . 

11.  Similar  to  figure  !0  but  with  .a  PC  disc;  cxn'pl  that  impact  h.i.s  .tlrcady 

occurred  in  frame  ,i.  The  PC  fails  catastrophically .  The  first  of  several 

ignition  si  te.s  can  be  detpcfed  in  frame  J.  Fic'ld  of  vi'.w  2.0  win.  Frame 

t  ir,ic7  tis  ;a,0  ;b  ,  34  ;  c  ,94  ;d  ,  I .'  I  ;t  ,  1 15  ;  f  , !  4 7  . 


12.  Experimental  arrangement  for  obtaining  pressure-time  curves.  W,  drop- 
weight;  H,  R!,R2,a3,  hard  steel  rollers;  C,  cylindrical  guiding  sleeve 
S,  sample. 

13.  Stress-strain  curves  for  samples  of  (a)  PP,  (b)  PC.  Catastrophic 
failure  of  the  PC  sample  occurs  at  a  strain  of  ca.  1.8.  No  discontin¬ 
uities  are  observed  in  the  PP  curve. 

14.  Experimental  arrangement  for  friction  experiments,  D,  detector; 

S,  slits;  CD,  chopper  disc;  SD,  sapphire  disc;  PP  polwer  pin; 

M,  balance  arm;  W,  counterweight;  DS,  drive  shaft. 

15.  Calibration  curves  for  the  friction  apparatus  (a)  high  temperature 
region  (b)  low  temperature  region. 

16.  Examples  of  the  traces  obtained  with  the  apparatus  (a)  is  a  soldering 
iron  at  5/5  g  (b)  a  small  diameter  PC  pin  and  (c)  a  large  diameter 
PC  pin.  H  is  the  signal  strength  when  there  is  no  filter  between 
source  and  detector  and  H'  the  strength  when  a  PMMA  filter  is  inter¬ 
posed  between  source  and  detector. 

17.  Polymer  pins  after  friction  experiments,  (a)  a  PC  pin  showing  many 
layers  x»’hich  have  been  sheared  off  and  built  up  on  the  trailing  edge 
(b)  a  PP  pin  which  has  rounded  and  shows  evidence  of  melting  (c)  a 
PMMA  pin  which  has  a  roughened  surface  caused  by  small  particles 
being  brittly  chipped  out  of  the  surface, 

18.  Micrographs  of  polymer  pin  surfaces  after  friction  experiments, 

(a)  PC  and  shows  considerably  strain  in  the  frictional,  direction. 

(b)  PMMA  which  has  a  powdery  appearance  with  many  small  particles 
chipped  out  of  rhe  surface. 
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The  temaihing  three  sections  of  this  renort  (3,4,5)  are  on  reaction 
kinetics  and  their  analysis.  A.11  have  been  published  in  the  open  literature. 

For  completeness,  reference  is  given  to  a  paper  (Chaudhri  et  al,  1977) 
based  on  work  performed  by  Dr.  Chaudhri  when  he  visited  Picatinny,  and  to 
other  papers  published  by  the  gr  'ncs 

M.M.  Chaudhri  and  J.E.  Field,  “Fast  Decomposition  in  the  Inorganic  Azides", 
from  Energetic  Materials,  Vol.  1,  F.d.  H.D.  Fair  and  R.F.  Walker  (Plenum, 

1977). 
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NEW  METHOD  TOR  TG  AND  DSC  DATA  ANALYSIS 
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All.NI  HAU  l 

The  paper  describes  a  graphical  computer  method  for  analysing  TG  and  DSC 
traces  v.liich  i'i\cs  ail  llircc  reaction  parameters  (A,  A,  ;ind  /;)  c!i;iracie!i>ii)g  an  «th- 
ordcr  reaction  from  a  single  trace.  If  the  reaction  proceeds  in  multiple  stages  (A,  A 
and  »)  triplets  c.in  be  obtained  for  each  stage.  The  computer  programmes  arc  basically 
simple  and  use  little  compt.iing  time  (typically  a  few  seconds).  The  advantages  of  the 
approach  over  earlier  methods  are  discussed.  As  a  test  of  the  method,  results  on  the 
dehydration  of  calcium  oxalate  monohydrate  arc  described. 

tNTRODl’CnON 

The  deduction  of  the  rca>..ion  parameters  E  (activation  energy),  A  (frequency 
factor),  and  n  (reaction  order)  from  TG.  and  DSC  traces  is  clearly  of  interest.  A 
number  of  methods  have  been  used  in  the  past  to  do  this  but  all  h,ave  unsatisfactory 
features.  In  the  work  described  here  a  number  of  Fortran  programmes  were  written 
to  extract  all  three  reaction  parameters  of  an  nth-order  reaction  from  a  single  TG,  or 
DSC  trace.  The  programmes  were  run  on  the  Cambridge  IBM  370, 165  computer, 
and  the  plotter  was  used  for  graphical  output.  So  that  the  advantages  of  the  present 
approach  can  be  appreciated  earlier  methods  are  brierly  discu.ssed. 

EARLIIR  Mt  ritODS 

Murray  and  White'  have  .shown  that  the  temperature  at  the  reaction  peak  of 
DTA  traces  can  be  used  to  obtain  E  and  A  from  a  number  of  traces  if  a  first  order 
reaction  (cqn.  1 )  is  assumed,  where  w  is  the  fractional  rcsidu.al  weight,  /  tlie  time,  T 
the  temperature  and  R  the  gas  constant. 


If  cqn.  (1)  is  differentiated  to  obiaiti  the  maximum  for  du  /dr  then  cqn.  (2)  can 
be  derived,  where  <l>  =  dl/df  is  the  heating  rate  and  the  temperature  at  the  rate 
maximum. 


In 


E  }_ 
R  ' 


(2) 


A  plot  of  )ir„  versus  In  yields  E'R  as  the  slope  and  In  (AR'E\  as  the  intercept 

of  the  straight  line  with  the  ordinate.  The  main  disadvantage  of  this  method  is  that 
it  only  uses  a  single  point  (r„)  from  a  trace  to  produce  one  point  on  the  activation 
energy  plot.  A  number  of  experiments  are  needed  to  obtain  a  reliable  value  for  £. 
There  is  also  the  limitation  that  the  theory  only  works  for  a  fust-order  reaction. 
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Kiisingcr*  mioaciI  Uial  the  icchnio/ic  «an  be  cxten.lcd  to  aa  /itii-order  rcaciiun 
by  deriving  cqn.  (3). 


=  A  n'C’ 


(3) 


He  shows  in  his  paper  that  the  prothict  n  tejl," '  is  independent  of  the  I.’cating 
rate  (0)  and  almost  equal  to  itnily.  J  herefore,  cqn.  (3)  cft'cciivcly  reduces  to  eqn.  ('2). 
The  approach  again  niahcs  poor  use  of  the  inforniaiion  containctl  in  a  IXFA  trace. 
Oraw  a's  analysis^  of  TG  traces  begins  with  a  general  form  of  reaction  equation. 


dw 

W 


(4) 


This  is  integrated  by  Ozawa  who  shows  that  the  rohition  to  the  right-hand  side  of  the 
equation  is  a  polynoniial  in  the  variable  JI/Rf.  Etiuation  (5)  is  an  approximate 
formula  which  relates  the  heating  rate  0  with  A*. 


log  0,  -b  0.4567  =  tog  0,  +  0.4S67  -  (5) 

where  T,  and  Tj^conespond  to  the  same  fractional  decomposition  («’)  of  the  sample. 
Since  a  plot  of  ibg  0  verstis  l/T  can  be  obtained  for  each  value  of  w,  Ozawa  aggregates 
all  these  results  in  a  master-curve  with  a  gain  in  accuracy  of  the  final  result.  Equation 
(5)  only  assumes  that  the  function  f(tv)  is  consliini  for  a  fciven  value  of  iv. 

Tliis  method  also  needs  a  number  of  experiments  at  different  heating  rales  to 
generate  enough  data  for  a  sin.dc  activation  energy  plot.  Although  the  method  was 
developed  for  ihcrniograviinciric  -ata,  it  can  be  applied  to  OSC  data  by  integrating 
the  peaks  to  give  values  of  tc.  A  theoretical  analysis  of  sucli  a  process  is  given  by  Reed 
ei  al.^. 

Borchardt  and  Daniels*  describe  a  method  which  allows  E,  A  and  n  to  be 
determined  from  the  sliape  of  a  DTA  curve.  Their  appi  oach  was  applied  to  reactions 
in  the  liquid  pliase.  The  type  of  reaction  to  whicl)  this  method  is  applicable  must  have 
a  single  rate  constant  and  the  activation  energy  must  not  vary  with  temperature. 
By  considering  the  heat  transfer  equation  of  the  c.ilorimcter  cqn.  f6)  was  obtained 


dt 


KA' 


"  dt  j 


where  U  is  the  number  of  motes  presem,  ,V„  the  initial  number,  K  the  heat  transfer 
cocflicicnt,  A  the  area  under  ihc  DTA  curve,  Cp  the  heat  capacities  of  the  two  liquids 
and  dr  the  icm|vrature  difference  between  the  two  cells.  Integration  of  cqn.  (6)  with 
respect  to  lime  gives 


=  {C,d7N.  An} 

where  o  is  the  area  under  tlic  DTA  curve  up  to  the  present  temperature. 
The  expression  for  the  rate  constant  A-  of  a  reaction  of  order  «  is 

■ 

where  V  is  the  volume.  Substituting  cqns.  (f.)  .and  (7)  for  dWIdt  and  gives 
r  C  ICAT 

k l" '  ”  ic  1 

1  No  J  {K^A'-a)  -c’Jry 


(7) 


(8) 


(9) 
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In  order  to  obl.iin  reaction  parameters,  Morcliardl  and  Daniels''  assumed  a  value 
for  H.  Tliis  allowed  k  to  lie  taleulalcd.  A  plot  of  In  k  versus  l/7‘  (activation  energy 
plot)  yields  a  straight  line  if,  and  only  if.  tlic  eoireci  value  of  ii  has  been  assumed, 
thus  values  of  /7,  A  and  n  can  be  obtained.  The  main  dilliculiy  of  this  method  is  the 
number  of  plots  that  have  to  he  made  for  each  value  of  ii  until  the  points  lie  on  a 
straight  line. 

A  number  of  computer  programmes  arc  presented  by  Benin  ct  al.‘  and  u.scd  to 
solve  three  problems,  namely  (i)  the  production  of  ihcrmograiViS  given  the  reaction 
parameters  E,  A,  n  and  given  that  the  reaction  is  of  the  H(h-ordcr  tyjie;  (ii)  the  calcula¬ 
tion  of  the  (E,  A,  ii)  triplet  from  an  cv|icrimcntal  thermogram  -  -  this  is  the  inverse  to 
problem  (i);  and  (iii)  the  determination  of  the  reaction  p.iinmeicrs  and  the  complex 
reaction  mechanism  from  a  given  c.xixrimcot.al  thcrmogr.im.  The  first  problem  is  a 
straightforward  function  plot  of  the  solution  to  the  nth-order  reaction  equation. 
The  second  problem  is  more  diilicnlt,  and  presents  itself  to  anyone  who  wants  to 
obtain  reaction  p.iramcters  from  experimental  data.  Benin  ct  al.*  treat  it  as  a  problem 
in  curve  fitting.  They  minimise  the  discrepancy  bcusccn  the  cxpcrimcnial  data  and  the 
calculated  function  by  taking  as  the  discrepancy  criterion  the  maximum  deviation  of 
the  function  from  the  data.  The  solution  of  the  third  problem  is  an  attempt  to  reveal 
the  mechanism  of  a  chemical  reaction.  Tlic  programmes  become  quite  involved  as 
they  work  out  activation  energies  for  each  section  of  the  thermogram  by  using  a  plot 
of  reaction  rale  versus  IfT  for  a  constant  amount  of  conycision,  a,  but  different 
heating  rates. 

The  most  genera!  form  of  the  fit  function  assumes  a  decomposition  process 
which  proceeds  in  a  numbei  of  individual  st.iges.  r.ach  stage  consists  of  an  elementary 
nth-order  reaction,  i.c.,cqn.  (1)  with  A,  ^and  ir  replaced  by  A^.  £(  and  w.-,  respectively, 
where  denotes  the  variable  belonging  to  the  /Ih  stage  of  tlic  reaction.  The  program¬ 
me  which  optimises  the  fit  of  this  function  to  the  dam  lakes  up  to  4  h  of  computing  time. 
This  makes  the  method  rather  expensive.  Another  weak  point  is  the  discrepancy 
criteria  between  the  data  and  the  fit,  which  is  taken  as  the  maximum  difference.  It  is 
well  known  that  the  data  points  towards  the  end  of  a  decomposition  process  become 
less  reliable  due  to  inipuiiiies  and  zero  errors.  But  these  points  have  the  same  weight 
in  the  analysis  as  all  the  others. 

Schempf  ct  al.’  describe  a  prognimmc  which  fits  a  least  square  polynomial  to 
the  weight-loss  curve.  The  quantity  dw/dr  s  calculated  by  diftcrentiation  of  the 
polynomial.  A  straight  line  is  then  fiitcvl  to  a  plot  of  log  k  versus  l/T,  and  the  least 
square  criterion  is  used  .again.  The  re.actioi  order  is  assumed  to  be  one.  and  no 
attempt  is  made  to  eliminate  the  deviations  l  .om  an  ideal  reaction  at  the  licgtnning 
and  end  of  the  weight-loss  curve.  Our  work  sf  ou  s  that  this  is  very  critical  and  c.in  be 
obtained  from  an  .Arrhenius  plot  performed  ly  the  eompiilcr. 

KEW  METHOD 

Tlic  method  developed  was  based  on  the  idc:  s  of  Borchardt  and  Daniels'.  Our 
data  were  obtained  from  both  a  ihermobalnnce  (St  mton  ReJerofe)  and  a  differential 
thermal  calorimeter  (Pei  kin- Elmer  DSC-2). 

The  analysis  assumes  a  rc.iction  which  '  m  be  described  by 


Since  the  heating  rate  h  (deg  min" '  is  constant,  the  time  variable  can  be 
eliminated  fioni  cqn.  (10)  by  putting  d7  'd/  -  h/OO  (deg  sec' '). 
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Taking  natural  logariihtns  on  both  sitlcs  gives 


,  .  ,  /I,  X  60 

In  iXiv)  In  — y- 


Note  tb.ll  ( -•  dii’/dr)  is  ;i  positive  qu.inlity  since  the  fr.iction.il  rcsulua!  weight  decreases 
with  inci  casing  temperature. 

A  plot  of  {ln(—  die,/.l/)  —  In  f.fn  )}  versus  {I/T}  only  yields  a  straielii  line  if 
f,{n’)  has  been  clioscn  eoireeliy.  In  the  rase  of /,(n')  —  m"  the  reaction  order  n  is 
varied  iin'il  .i  straight  lino  is  found.  The  activation  energy  /'can  be  obtained  from  the 
slope  of  this  plot,  the  fre.gioncy  factor  ,1  from  the  intercept  on  the  ordinate  and  (be 
reaction  order  //  is  determined  by  the  strai.shl  line  critciion. 

This  analysis  can  .ilso  be  applied  to  the  D.SC  data  since  the  basic  r-qn.  (12)  only 
differs  from  (10)  by  constants. 


=  ±  -  II-,)  .4  exp  ^ f(«') 


The  i  (  — )  in  the  equation  corresponds  to  an  cxoihciinic  (endotherinicl 
reaction  where,  ()  is  the  he.il  given  ofl'by  the  sample,  d//  is  the  cnih.ilpy  change  of  the 

entire  reaction,  u’o  is  the  initial  weight  and  >r,  is  the  re:.l  weight,  and  the  other  para¬ 
meters  arc  as  defined  previously. 

Taking  logarithms  gives 

In  j  ~  •*'  “  «■.)  ^  ^ 

This  time  dQ/iU  at  temperature  Tis  the  mcajnred  quantity  and  «■  needs  to  be 
calculated.  This  is  easily  done  since,  if  the  reaction  results  in  a  weight  loss 


‘  f  ^ 


The  computer  numerically  integrates  dOfdt  to  yield  the  weight  at  lime  {.  A  plot 
or(lndg/dr  -  hi  f(»)}  versus  1/T produces  r.  d  and;;  (iff(i.  )  =  w")  in  the  same  way 
as  for  the  wvight-loss  curves  above.  I  he  slope  is  again  - /I/R  and  the  intercept  is 

lnd//(ivg  - 

In  practice  tlic  data  from  the  TG  or  DSC  is  recorded  or.  8  track  ASCI  paper 
tape.  This  inform.ilion  is  then  read  into  disc  files  of  the  computer.  These  files  can  be 
used  by  the  programmes  which  c.ilciilaic  the  acrivalion  energy,  etc.  .^flcr  these 
programmes  have  run,  the  plotter  output  is  viewed  on  a  ivIev  ision  screen.  Only  if  the 
plot  is  salisfaciory  is  a  copy  obtained  from  the  plotter  for  the  final  analysis.  The 
dimensions  correspond  to  the  origiiwl  o-jiput  of  the  chart  recorder  so  that  a  direct 
comparison  is  possible.  DilTerent  pen  colours  arc  u-ed  to  dincrcntiaie  between  the 
raw  data  (black)  and  the  ilicorctic.il  fit  to  the  dat;  (green),  these  colours  are  re¬ 
presented  in  Fig.  1  as  .1  broken  line  (lilacK)  and  n  full  line  (green). 
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I  in.  1.  Wii!;lil  ’.O';',  tvirvc.'  of  Cci  C"0!lli’f)  HnaK^is.  The  d.isiK’t)  )iiie  r.  the  oiieiiv.!  iliit:ian<l  Ihcfiiil 
line  ihc  lei'h't  iisiri,".  ;l>c  r.':.  pnicJ  if.,  .i,  ■/')  iiijilci  fc-i  the  i!»ce  ■>iin!;'.iii  Inx  ee'Ctioos  illicit  arc 
s<.-l>iiia;cv'  I';.'  vouic.'ti  ■.Kl^!\e^!  iiiics. 


The  {)i/fci  cni!al  pioi^irinwic* 

'I  itcrc  urc  (htcc  ctilTcrcnl  'v.iys  in  v-Viiclt  tiiircvcniinl''  to  i!\c  v.ci.ftl’t  los'i  ciir\c 
were  obUiiiiCil.  The  firsl  uses  ili!tVri:n.:cs  niul  i>ri>iluccs  a  noisy  tiifTciciilial  eui  vc.  xfliich 
is,  however,  nenirsl  to  ihc  criyin.il  d.-tl.-!.  The  other  two  u.ve  a  polynomial  lit,  and  a 
sirliiic  (it,  respect ieely.  The  poiynotnial  f.i,  n.'.'.ally  ol  25ih  cndcr,  is  yeneiMlIy  yood  bitl 
prodnec.s  artilicially  laryc  innini.itnms  over  sivaieht  regions  of  the  curve.  The  best 
reviilts  were  ohlaincd  by  a  citioc  sphne  hi.  This  jsioyrarnine  fits  a  third  oidcr  curve  to 
sections  of  the  data,  and  joins  ilicni  in  sti-.h  v,a>  titai  titc  function  and  Ihc  first  two 
derivatives  arc  continiu-ns.  lii-  procr  i-.n.nc  iis;s  a  si.tndard  siibrouline  from  the 
Harwell  !  oil  ran  Mtraiy,  e.  in,  h  ,  lntiv-.es  !i\  ouii !  i:..;s  (points  where  two  sphnes  meet) 
accordiiif.  to  the  b.  havi'nr'  of  '  .-.riein  ei' ve  i  iiis  means  that  over  a  slraiplit  part 
of  the  Cifnc  a  lo.np  section  '  .'in  !>e  litic'i  ov  a  .i  t;  le  euhiv'  spline,  wljcrcti.s  at  places  of 
greater  variation,  smaller  setiions  an  ebi’Sct  to  Itc  appro.vimaled.  The  diflVrential 
of  the  curve  IS  i’.cn  ..Ii’-tined  by  lire  tltird  order  spline.  The  fit  of  tlie 

spline  and  i!’c  uifTeiciniai  is  wrv  aci;n;a'.c  The  tmomit  of  data  vmoolhina  nclhevcd 
by  this  r.i  suits  the  purpose  oi  iiio  folieo' iae  paiera/nine  (iic;>,t  section)  best. 

The  f>n>p/:iinr,h-  w  ftn.f  F'.  i  am:  n 

It  is  assumed  in  this  an.'ilvsi.s  that  the  iT:iriion  that  leads  to  the  weight-loss  curve 
is  at  least  piccc-wtsc  of  the  nth  order,  i.n.  follows  eqn.  (10)  with  (((;>')  ~  it"s  This 
prograniine  plots  [In  (—  du  'd.')  -  n  In  re;-  versus  itT  for  dilTercnt  values  of  n.  The 
values  of  dic/Jt  and  of  n  ;;re  obtained  from  a  tile  that  was  created  by  the  previous 
prof.iamtnc  .and  <ont:iins  the  .smoothed  wciehi-losv  data  and  tire  difl'erenlial  in  the 
form  of  Ihc  .spline  parameters.  The  value  of  ii  is  taken  as  zero  for  the  first  plot  and 
liicn  incteasci.l  i-.y  (1  ."I  ren  times  up  to  a  ■-  2.  I'hesc eleven  lines  h;ii c  different  cui'vaturcs. 
The  straii’htcst  of  them  tan  easily  be  selected  by  eve  ami  tins  n  vr.ltic  is  then  taken  as 
the  reaction  ore.ei  Fh-.’  -iope  ytiMs  /  .md  tie.  uucucpi  ri.  ihc  rc.suitin''  triplet  of 
mimbcrs  (AJ,  .(,  /;)  ntc  newessary  aiid  -.uricieni  to  eliaiacterisc  an  mh-ordcr  reaction. 
The  type  of  icactir’ii.  i  e  f(ic)  -  u*  e.m  be  chaueci.t  c.isdy  so  ihai  ilx  puipiainnic  is 
apivheablc  to  reactions  of  any  tv  pc.  i.e.  f(e  )  tan  have  any  form  (c.k.  ramlom  chain 
scission). 

This  and  d;e  fuiti.mi.ii'  |M.-y.iai’)ir,i-  rrm  I-..-  ot-l  iiiml  fu.iv,  ihi  rmlhors. 
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The  pro^'jammc  to  pn.'Jitce  a  Ihcorcfiral  fif  funJt  ftu*  ( A,  .-t,  tt)  trtplet 

Tlii'i  profrainnie  allows  the  production  of  a  wci;i(it-loss  ciirrc  niih  up  <o  four 
diffcrciil  sections  e;u!i  of  « liich  has  iis  own  H.  A  and  n  factor,  ll  was  iutpl.-iuentcd  on 
a  Hcvvliit  Packard  1  >  :sk  lopcalcidalor  and  ploiir  r  of  the  l>pe  HP  V8;J5.  The  ihcoiciical 
plots  produced  in  such  a  way  fit  the  ori-inal  wciflu-loss  cui\e  very  well,  fhey  al>(» 
provide  a  useful  cherk  that  at!  the  cakufaiio.is  necessary  to  pioducc  this  triplet  were 
correct,  and  pioducc  a  pood  fit. 
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t  ig.  2.  Atiiv.aliou  cncipy  (ilol  for  vsiiotis  values  of  roaeiinn  oidt-i  ».  The  n  t.-ihic  which  f.itcs  IIk^ 
bcsl  .sliaixiil  line  M':’.i''..'n(s  is  the  one  iipproiviale  (o  the  u.ieiion  (here  0  6)  liiere  ere  Ihiie  sups  lo 
the  reaction  but  ibe  lust  two  coricsiiond  to  only  ca.  2,0  •)(  ihc  weight. icis. 


II  ST  OF  MI.T1101J 

The  dcliydralion  reaclioir  of  calcium  oxalate  inonohydrate 
Ca  CjOaHjO  Ca  CjO*  i  H,0 

is  well  understood  and  values  for  the  reaction  p.iramciers  can  lie  obtained  from  llie 
literature*'  It  therefore  provides  a  good  test  c.isc  fisr  tl-.c  new  analysis  method. 
Figure  I  shows  (he  wcichi-loss  ciii  se  of  a  saiupie  of  Ca  C. 04(1,0. 1  ho  dashed  curve 
represents  the  original  dat.i  and  the  ftil!  ci  rvc  the  theoretical  fit  based  on  the  computed 
(C,  A,  n)  tiiplcl.  Figure  2  shows  titc  a-,  (iva'ion  cuercy  plot  The  number  at  Ihc  cn<* 
each  line  gives  the  rc.'iction  order  ih.al  is  assunird  fi'r  the  calculaiion  of  the  cr'iic 
spending  line.  The  'straight  lines  which  nppioxim.ric  the  cui  ves  best  arc  used  to  obtain 
(£,  /I,  n)  tripl.'ls  for  the  ;hrcc  sections.  Since  these  (/T.  .1,  w)  triplets  completely 
Specify  the  reaction  the  origiral  weightdoss curve  c.in  be  modelled  (full  curve  Fig.  1  k 
The  results  obtaincvl  with  this  method,  and  iitcrature  values  arc  sunimari/ed  in  Table 
I.  All  our  esperiments were  performed  at  a  heating  rate  of  10  dec  min' '.  The  results 
obtained  in  this  way  have  an  error  of  -  5  kJ  mol  ‘  ’  w  hich  is  doc  to  the  variance  of  the 
samples  themselves.  The  method  itself  has  much  gtcaicr  as'curacy.  It  was  found  that 
the  best  straight  line  could  be  fitted  to  a  curve  corsesponding  to  a  reaction  order  of  0  6 

TABLE  t 
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(kJ  mole-') 
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which  indk-atcs  a  surface  rcaclioii  (Iheorclically  2,0).  Il  can  aKo  he  .seen  on  Fig.  2 
that  tiicie  are  three  (.lislincf  stages  of  the  ivaelion.  T))c  lirst  iwo,  hoaeecr,  corre>ix<nd 
to  only  2%  of  the  eight  and  ate  attributed  to  vtthiile’  substances  on  the  surface. 
This  s)k'«  s  the  scnsiiixity  of  the  method  used  cs|Kaially  considering  that  the  \crtical 
axis  has  a  logaritlimic  scale.  The  teal  dehydration  reaction  is  rcpivscnted  by  the  Iasi 
straight  line  .segment. 

The  method  outlined  in  this  paper  has  recently  been  successfully  applied'®  to 
decomposition  studies  of  f’.F.T.N,,  P.F.T.N.  with  fillers,  liigh  densily  polyethylene 
(IIDPE)  anti  polytctrafluorocthylcnc  (I’TFL). 
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SECTION  4 


ANAl.YSIS  OF  ISOTHERMAL  KINI  TK'  DATA 
FROM  SOLID  STAI  C  RFACTIONS 

Tonc  B.  Tang  and  M.  M.  Ciuijdiiri 

Phyiles  and  Chemistry  of  Solids,  Cut'rmlish  Lahoralory,  Madinghy  Road,  Caihbtidye  UK 
(Received  Dceeml'ft  12.  1978) 


In  mosi  solid  slate  icactions  the  reaction  vcloiily  van  be  described  as  a  product  of 
two  fnnciions  A(7  )and  /<!  -  j)  wheic  7  is  ihe  temper.iinrc  and  <  the  siegret  of  con¬ 
version  of  the  solid  reactani.  I  he  physical  interpretation  of  tWese  functions  is  discussed, 
and  a  systematic  metnod  is  dcscritM.’d  by  winch  ft  1  j)  of  a  reaction  i«  identified  from 
its  kinetic  data.  A(T)  and  the  icaction  metbamsm  arc  then  determined.  1  hi  '  od 
has  been  successfully  applied  to  analyse  the  kiiretic'  of  the  iheinul  itcconipo  .n.  n  of 
silver  azide. 

In  a  solid-state  rcactinn,  the  re.'tciiyyn  velcrity  is  eivcn  by  -  d/dt  (1  —  «)  »  a, 
where  a  =  a(f  )  is  tlie  fraction  of  the  solid  reaclant  wlntli  has  reacted  by  lime  /. 
Its  kinetics  can  be  solved  by  dctcrminiii);  i  .as  a  function  of  iem)veraturc  and  the 
global  amount  of  reactant  left.  This  phenomcnoK'uical  knowlcdtir  is  n  necessary, 
though  not  sullicient.  condition  f<»i  elucidating  the  reaction  mechanism.  U  is 
necessary  to  formulate  the  reaction  velocity  in  terins  of  the  global  variable  a, 
because  there  is  a  continuous  collapse  of  structure  in  tiic  reactant  Furthermore, 
the  local  concentration  of  ihc  reactant  varies  throughout  the  reaction  volume  and 
cannot  be  used  as  a  state  variable.  In  fact,  unlike  the  case  of  a  homogeneous  re¬ 
action  in  the  liquid  or  gaseous  phase,  there  is  no  real  'reaction  order'  with  respect 
to  any  reactant  in  a  reaction  involving  condensed  matter  whose  mechanism  is 
usually  of  the  heterogeneous  type. 

If  the  reaction  proceeds  isothcrmally,  it  is  observed  empirically  that  (a.  i) 
curves  corresponding  to  different  temperature.v  T  arc  isomorphic  to  one  another, 
at  least  within  a  range  of  T,  j.e.,  by  a  linear  scale  change  in  r,  dilTcrcnl  curves  can  be 
superimposed  [1],  It  follows  that  »  is  a  separable  function: 

(*)i«>lhtfn:«l  "  ft  I  ~  •  (1} 

Hcrc/^l  -  a)  may  change  in  dilTerent  ranges  of  T  or  a.  For  every  /(I  -  a),  there 
corresponds  a  single  A'(7').  It  should  be  noted  that  experimental  data  may  be  ade¬ 
quately  analysed  by  U)  only  if  it  has  been  ensured  that  tlie  tcmpei  aturc  distribution 
in  the  sample  is  su/licicnily  uniform  and  corstant.  rurihcrmorc,  the  theoretical 
significance  of  A(r)and  /i;i  -  aldcterrnined  froi  i  the  data  should  always  be  c.xam- 
ined  wiihj  icgard  to^  the  class  of  mechanisms  .hey  indicate.  A  question  of  con- 
tstency  arises  in  this  respect.  In  the  literature,  coherent  and  integrated  .accounts 
of  the  physical  fas  contrasted  with  the  formal)  meanings  of  A'(7')  and  /{I  -  a) 
are  not  easily  found.  A  discussion  on  their  ir.erprciatic'nv  therefore  forms  Ihc 
first  part  of  this  paper.  In  the  second  part  we  propose  an  ellicicni  meihvxl  of  de¬ 
termining.  with  confidence,  both  A I  -  aland  A  7’ l  from  y  and  »(/)data.  The  inter¬ 
pretation  of /(I  —  »)  is  the  essential  link  in  this  method  between  the  experimental 
data  and  the  functions. 


Physical  inicrpreta  ions 


The  fiwriioii  f(l  -y) 


Solid-state  rc.aetions  are  complex  pri’ccssv  -  which  proceed  in  scvcial  si.iges. 
These  can  be  the  dvlocaluation  or  ir.in  .lci  o!  .  lections  in  cl.cmical  bonds  (in  the 
ease  of  lum-mctals),  the  dilfiision  of  aioms,  lie  i.idicals,  or  ions,  the  desorption 
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of  product  molecules  when  they  are  in  the  gaseous  phase,  the  heat  transfer  to  the 
reaction  zone  in  the  case  of  endothermic  reactions,  and  the  formation  of  a  new  solid 
structure  (crystalline  or  amorphous)  if  one  of  the  products  is  in  the  solid  phase. 
The  last  step  may  often  be  further  differentiated  into  nucleation.  growth  of  nuclei 
(at  velocities  which  depend  on  sizes  of  nuclei)  [1]  and  sometimes  the  collapse  of 
the  lattice  from  a  transitory  one  to  the  equilibrium  structure  [2].  Irrespective 
of  the  details  of  reaction  mechanisms,  however,  under  a  given  set  of  circumstances 
(T.  x.  sample  history  etc.)  one  of  the  stages  will  be  the  slowest.  It  then  acts  as  the 
rate  limiting  step  of  the  reaction,  and  it  will  determine  the  kinetics  i.e..  the  rate 
lavs  in  ( I ). 

Further,  a  solid-state  reaction  has.  in  contrast  to  a  homogeneous  reaction  whose 
progress  is  independent  of  spatial  coordinates,  an  additional  controlling  factor, 
namely  topochemistry.  This  refers  to  the  geometrical  shape  of  the  solid  reactant 
and.  in  different  cases,  to  its  free  surface  area,  its  defect  structure,  the  thickness 
of  the  product  layer  if  solid,  or  to  the  product- reactant  boundary,  etc. 

The  function  1  —  a)  reflects  the  nature  of  the  rate-limiting  step  and  the  topo¬ 
chemistry  of  the  reaction.  It  may  accordingly  depend  on  certain  sample  conditions, 
such  as  whether  the  sample  has  been  pre-irradiated  or  bleached,  and  whether 
the  sample  is  in  the  form  of  a  powder,  or  a  large  single-crystal  of  a  different  shape 
from  the  crystallites.  It  will  vary  in  several  ranges  of  T  if  in  each  of  them  a  different 
elementary  step  becomes  rate-limiting,  as  occurs  in  the  decomposition  of  potassium 
azide  [3].  At  a  given  temperature,  it  may  also  change  in  different  ranges  of  a. 
due  to  the  switching  of  the  rate-limiting  step  or  topochemical  changes.  This  hap¬ 
pens,  for  instance,  in  the  oxidation  of  zirconium  [4]  and  in  most  decomposition 
processes  [1].  An  extreme  case  is  the  decomposition  of  ammonium  perchlorate 
which,  in  the  two  temperature  regimes  below  and  above  f>lQK,  has  entirely  dif¬ 
ferent  reaction  mechanisms  and  in  fact  yields  different  reaction  products  [S], 
In  such  cases  there  may  be  competing  paths  for  the  chemical  reaction  or  it  may 
in  fact  be  followed  by  another  chemical  reaction  whose  'onset  temperature’  is 
higher.  In  all  possibilities,  however, /T1  -  a)  should  be  the  same  for  a  given  a 

incependent  of  /'  (within  a  range)  if  it  is  to  have  more  than  only  an  empirical 
significance. 

In  Table  1  we  have  collected  together  the  more  common  forms  of /(I  -  a) 
which  have  been  used  in  the  literature,  and  the  corresponding  integrated  forms 

\  df  iif(\  -  x)  =  j  dx/yil  -  *)  =  f(»).  Note  that  /'(a)  ^  K  [i  ~  t„)  if  the  range 
0  0 

of  X  for  which  it  becomes  applicable  starts  at  x  *  x(fn).  Also,  for  simplicity  here- 
aAer  we  write  K  for  K{T). 

In  many  reactions,  such  as  most  dectn-f  csitirrs  and  dehydrations,  the  rate- 
limiting  step  takes  place  at  the  imerflce  between  different  phases  as  in  sublimation. 
The  speed  at  which  the  interface  icbctant  is  (at  a  given  temperature) 

then  either  a  constant,  or  a  unique  function  of  the  interfacial  area.  This  area 
therefore,  from  the  kinetic  point  of  view,  plays  the  same  role  as  that  of  con¬ 
centration  in  homogeneous  reactions.  If  the  speed  is  constant,  then  the  theoretical 
significance  of /(I  -  x)  is  clear;  it  gives  the  area  expressed  as  a  fraction  of  the 
original  area  at  x  »  0.  This  is  the  case  of  a  reaction  controlled  by  the  movement 
of  a  coherent  phase-boundary  and  listed  as  C  and  H  in  Table  1.  In  this  situation, 
the  explicit  form  of/ll  -  x)  depends  on  the  geometry  of  the  reacting  system, 
though  generally  it  is  a  decreasing  function  of  x  or  at  most  constant. 

If  the  reaction  consists  of  the  formation  of  compact  nuclei  of  a  solid  product 
at  localized  places  in  the  reactant  followed  by  their  relatively  rapid  growth,  then, 
to  express  the  total  interfacial  area./[l  -  a)  is  derived  from  the  laws  of  nucleation 
and  growth.  This  is  the  situation  when  the  reaction  is  autocatalytic  (13}:  reactant 
molecules  at  a  reactant-product  interface  react  in  preference  to  those  at  a  reactant- 
'  vacuum'  surface.  The  preference  is  due  to  the  existence  of  microstrains  in  the 
reactant  at  the  interface,  or  due  to  the  electrochemical  potential  of  the  product 
phase  when  the  rate-limiting  step  is  a  redox  process.  The  various  possible  forms  of 
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/( 1  —  i)  for  an  autocatalytic  reaction  are  listed  in  Table  I .  ^4  to  £.  Note  that  they 
have  the  general  form  /(I  -  s)  =  ^'“(1  -  a)".  In  .4.  B  and  £.  q  is  zero  and  a 
increases  monotonically  with  a.  Such  a  situation  is  most  unlikely  to  last  up  to 
a  =  I.  These  types  of  /(I  -  a)  therefore  may  apply  only  to  the  acccleratory  part, 
if  it  is  present  in  the  a-time  curve  and  which  will  usually  be  followed  by  a  decay 
part.  In  C  and  D.  </  is  non-zero  and  these  types  give  sigmoid-shaped  curves.  The 
inflexion  point  occurs  at  a'  =  pi{p  -f  q).  as  can  be  seen  at  once  from  the  condition 
a  =  0.  It  is  thus  kinetically  feasible  for  them  to  fit  the  complete  e.xperimental  curve. 
There  are.  however,  physical  grounds  to  consider  that  even  they  should  be  used 
10  analyze  only  the  acceleratory-  period  [1]. 

In  other  reactions  the  rate-limiting  step  is  not  confined  to.  or  does  not  only 
occur  ai.  the  reactant  surface.  For  instance,  in  the  unimolccular-decay  type  of 
reaction,  all  molecules  whether  on  the  surface  or  in  the  bulk  have  an  equal  proba¬ 
bility  per  unit  time  of  reacting.  This  is  the  ca.se  when  the  change  from  the  re¬ 
actant  to  the  solid  product  phase  involves  little  re-arrangement  of  the  reactant 
atoms.  The  reaction  has  homergeneous  mechanism  and  thus  a  true  reaction 
order  of  one.  i.e.  /( I  -  a)  =  1  —  Many  decomposition  reactions  tend  to  this 

Table  1 

The  common  types  of  solid-state  reaction 

Reaction 


Autocatalytic  Power-law  nucleation  and  growth  at  constant  speed:  A 

Linear  branching  chain  of  nuclei,  no  overlap  during  growth:  B 
Branching<hain  nucleation.  interference  during  growth:  C 
Random  nucleation,  growth  accompanied  by  ingestion  of  nuclei:  D 
Instantaneous  nucleation.  size-dependent  growth:  E 

Phase-boundary  controlled  decay  I -dimensional:  F 

2- dimensional:  C 

3- dimensional.  H 

(Jnimolecular  decay :  I 

1 -dimensional:  J 

Diffusion-controlled  2-dimensional:  K 

3-dimensional:  L 
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i)  Log  [  —  In  (I  —  »))  =!  Constant  +  m  log »,  0.15  <  i  <  0.5 

ii)  Plot  of  L.H.S.  in  i)  against  lug  t  distinctively  concave  upwards:  C,  and  arc  constants 
while  CfD  is  function  of  temperature 

iii)  h  =  0.774.  0.700.  0.664,  0.642 .  0.556  for  n  =  2.  3.  4.  5 _ x 

is  )  Alternative  derivations;  wc  have  obtained  the  approximate  forms  of  ft  I  —  *)  by  ex¬ 
panding  into  scries  (I  —  and  (1  -  »)■•''  to  second  order  in  i:  resulting  error  as  6 
t<10*;  for  7  <  O.S) 


37. 


limit  at  high  a  values.  Another  example  is  when  the  rate-limiting  step  is  the 
migration  of  product  ions  along  the  dislocation  network  to  form  additional 
growth  nuclei  at  dislocation  nodes  [17].  This  leads  to /[I  -  a)  =  x,  the  same  form 
as  for  branching  nuclei  [7].  The  slow  process  of ‘ageing’  in  some  explosives  when 
they  are  stored  at  room  temperature  may  be  by  such  a  nucleus-chain  mechanism. 
A  third  category  consists  of  reactions  controlled  by  the  diffusion  of  reactants 
across  a  product  layer  which  is  solid.  The  diffusion  may  proceed  uniformly  through 
the  bulk  of  the  layer  and  is  thus  structure-insensitive,  or  preferentially  along  its 
gross  lattice  imperfections  arising  from  product-reactant  mismatch.  In  the  case 
of  uniform  diffusion,  the  speed  at  which  the  product-reactant  interface  moves  is 
a  function  only  of  the  product  thickness  (and  temperature),  and  the  appropriate 
forms  of  /fl  —  a)  are  included  in  Table  1  as  7,  Af  and  L.  The  oxidation  of  metals 
often  follows  diffusion-controlled  kinetics:  in  sheet  form  these  tarnish  according 
to  the  parabolic  law  a'  k/.  Exceptions  aret  hose  metals  in  Groups  la  and  Ma  of 
the  Periodic  Table.  Excluding  beryllium,  they  all  form  oxide  layers  which  are 
porous,  so  that  the  atoms  of  the  metal  do  not  have  to  diffuse  through  a  coherent 
layer  before  coming  into  contact  with  oxygen. 


TJje  function  K(T) 

It  is  almost  always  the  case  that  the  temperature-dependent  part  of  (1)  can  be 
represented  successf^ully  by :  - 

K{T)=  K^exp(-Eikr)  (2) 

in  which  k  is  Boltzmann’s  constant,  and  the  macro-kinetic  constants  E  and  K-r, 
do  not  depend  on  T  (within  the  range),  though  usually  they  take  on  different  values 
when  /( 1  -  a)  changes. 

If  it  is  established  that  the  reaction  is  rate-limited  by  a  diffusion  or  migration 
process,  the  interpretation  of  KfT)  is  complicated,  but  obviously  it  is  proportional 
to  the  corresponding  transport’coefficient,  which  in  general  is  an  exponential  func¬ 
tion  of  r.  An  over-simplified  theory  for  the  situation  of  uniform  one-dimensional 
diffusion  gives  KiT)  =  (S/ Kn)*£)(T),  where  S  is  the  interfacial  area.  the  initial 
volume  of  the  reactant,  and  D(T)  the  diffusion  coefficient  (cf.  [19]). 

For  a  single-solid-reactant  reaction  which  is  controlled  by  a  surface  process, 
on  the  other  hand,  the  simple  theory  of  Shannon  [18]  is  often  successful.  This 
theory  is  a  generalization  of  the  Polanyi  —  Wigner  equation.  Assuming  the  existence 
of  some  activated  complex,  which  as  a  transition  state  can  be  treated  in  thermo¬ 
dynamic  equilibrium  with  the  reactant,  be  related  the  pre-exponential  factor 
to  the  rotational  and  other  internal  degrees  of  freedom  of  a  reactant  molecule 
in  addition  to  the  vibrational  ones.  Following  Shannon  we  can  set;  - 


A.  =  (A:77A)exp(dSVA) 

(3) 

and 

exp(  — £/k7’)  =  AH'jkT). 

(4) 

Here  the  mean-frequency  factor  kTIh  containing  Planck's  constant  is  usually  in 
the  region  of  lO**  s" '  (see  below).  dS*  and  dW*  are  respectively  the  entropy  and  the 
enthalpy  of  formation  of  the  transition  complex,  <5  is  the  thickness  of  one  mono- 
layer  and  V„  the  initial  volume  of  the  reactant,  and  S„((l  -  a)  gives  the  free  sur¬ 
face  or  the  product-reactant  interface  area  when  the  degree  of  conversion  is  a. 
(Strictly  speaking,  it  has  been  assumed  that  the  reaction  proceeds  isobarically). 

Note  that  in  this  interpretation  the  empirical  quantity  contains  the  surface- 
-to-volume  ratio  and  so  depends  on  the  sample  geometry.  Also,  it  is  apparently 
proportional  to  T.  (In  gas  reactions,  the  collision  theory  gives  x  T*.)  In  our 
opinion,  however,  if  the  vibrational  modes  are  being  considered  then  only  at  low 
temperatures  will  the  peak  distribution  of  phonon  frequencies  lie  at  kT/li.  For 
most  substances  (with  the  exceptions  of  Be.  Grand  diamond)  the  Debye  tempera¬ 
ture  is  less  than  500K.sothat  the  frequency  factor  should  stay  as  A  tip  /t  <  10'*s*' 
for  all  likely  experimental  temperatures. 
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The  factor  exp(4S'  Vt)  may  alternatively  be  written  in  terms  of  partition  func¬ 
tions  as  Q^!Q.  which  can  be  determined  from  spectroscopic  data  [20].  In  most 
cases  AS*  cannot  be  larger  than  the  reactant  entropy  of  melting,  and  exp{AS^Ik) 
comes  out  normally  between  unity  and  10*.  Occasionally  np{AS*ik)  is  found  to 
be  less  than  unity,  as  is  the  steric  factor  in  gas  reactions.  Such  a  negative  value  of 
AS*  means  that  the  activated  complex  is  more  ordered  than  the  reactant  (e.g. 
[21]).  Experiments  on  some  decomposition  reactions  have  given  K  which  are  a^ 
normally  high  in  comparison  to  the  theoretical  values  of  (3).  Hypotheses  put  for¬ 
ward  to  explain  such  discrepancies  include  co-operative  activation  [22],  proton- 
delocalization  [23].  and  a  mobile  layer  of  molecules  on  the  reactant  surface  [18]. 


Kinetic  analysis 

Current  practice 

In  determining  the  kinetics  one  wishes  to  find  K and /][  1  —  x)  or  equiv¬ 
alently  f(s)  so  that  the  reaction  velocity  can  be  predicted  at  any  given  T  and  a. 
This  is  commonly  done  by  analyzing  a  set  of  a(r )  or  equivalently  i(t)  values  ob¬ 
tained  by  monitoring  a  number  of  samples  reacting  isothermally  at  a  number  of 
temperatures.  The  consistency  of  the  and  £  values  with /HI  —  a)  should  as  far 
as  possible  be  assessed,  and  correlated  with,  for  instance,  microscopy. 

A  quick  method  of  calculating  £  was  used  by  Haynes  and  Young  [24].  Consider 
a  set  of  (a.  /)  curves  which  have  been  found  to  be  isomorphic.  For  any  two  curves 
(a,,  ^i)  and  (a.^,  4)  corresponding  to  temperatures  T,  and  respectively,  one  can 
write 

f(a,)  =  t,/Cr  exp  (-£/AT,)  ^5^ 

£(aj)  =  fj/C,  exp{-£kr.i). 

By  choosing  points  corresponding  to  the  same  a  on  the  two  curves  so  that  /l(a,)  = 
=  £(a2).  £  can  be  evaluated  by  plotting  In  /  vs.  1  T.  On  the  other  hand,  to  deter¬ 
mine  £(a)  often  a  trial-and-error  method  is  resorted  to.  Conflicting  forms  of  the 
function  have  sometimes  been  asserted  by  several  authors  for  the  same  material, 
like  NH4C10<  (see  [25])  and  KM„0,  (see  [26]). 

A  conventional  way  of  superimposing  isothermal  curves  is  to  convert  them  into 
‘reduced-time  plots'  by  individually  scaling  their  f-axis  with  the  factor  l/ti.  where 
Tj  is  the  time  when  a  =  0.5  on  the  /-ih  curve.  In  this  way  £(£()  is  absorbed  into 
each  scale  factor  and  all  (a,  t)  =  (0.5.  1)  points  coalesce,  while  other  a(r)  points 
may  be  plotted  out  to  see  if  the  curves  are  indeed  isomorphic.  Sharp  et  al.  [27] 
tabulate  the  theoretical  values  of  a  against  t'r  for  some  of  the  £(a)  shown  in  Table  1. 
They  propose  that  by  comparing  experimental  data  with  such  master  values  the 
correct  £(a)  can  be  identified. 

The  above  method  may,  however,  result  in  ambiguity  due  to  a  number  of  as¬ 
pects.  Experimental  data  contain  random  errors,  but  no  simple  statistical  analysis 
can  be  applied  to  the  identification  criterion  it  employs  because  no  straight-line 
graphs  are  involved.  Additionally,  a  general  problem  for  all  isothermal  experi¬ 
ments  is  the  zero-time  uncertainly.  The  finite  lime  taken  by  the  sample  to  reach 
the  designated  temperature  may  be  negligible  relative  to  r.  yet  may  affect  the  com¬ 
parison  with  the  tabulated  values  [28].  Moreover.  Ffx)  may  change  in  different 
regimes  of  the  <t(t)  curves,  as  mentioned  earlier. 


A  new  method 

Here  we  suggest  a  step-by-step  approach  to  determine  £|x).  It  was  noted  by 
Hancock  and  Sharp  [28]  that  for  many  forms  of  f(a|.  the  plot  of  log[-  ln(l  -  »)] 
Dj.  log  t  is  almost  linear  if  x  is  restricted  to  between  0.15  and  0.5.  Using  a  com¬ 
puter  program  to  generate  artificial  values  and  their  log-ln  plots,  we  have  found 
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this  true  for  all  the  theoretical  forms  listed  in  Table  1.  \sith  the  exception  of  B. 
C  and  E.  The  slope  in  each  case  is  listed  there  under  the  Column 'm  . 

Obviously  a  log-ln  plot  is  not  very  sensitive.  If  we  were  to  rel>  solely  on  it  to 
discriminate  between  the  functional  forms  of  Fin),  the  experimental  data  would 
have  to  be  of  the  highest  quality.  A  slightly  more  sensitive  way  is  to  plot  (  —  ln(  1  — 
-  il]'  IS.  I.  but  then  m  can  only  be  obtained  iteratively.  Fortunately,  m  does 
differ  significantly  between  different  groups  of  Ftx).  and  the  final  discrimination 
is  easiK  .ichiesed  by  a  further  graphical  step.  There  are  three  possible  situations 
tor  this  second  step;  - 

111  2  or  lo'.'-ln  concave  upwards 

It  will  be  seen  from  Table  1  that  this  situation  suggests  an  autocatalytic  reaction. 
U'r  which  j  Ai>''(l  -  11“'  for  certain  p  and  <y.  From  the  experimental  data  of  i 
.Old  /  one  can  then  do  a  least-squares  lit  on  the  graph  of  1  log  •'  (  log  y  against 
-J  log!  I  -  1)  .1  log  i.  and  find  />  from  the  r-intercept  and  q  from  the  slope.  Flere 
I  log  i  =  log  i(/|)  -  log  etc. 

For  an  i  expression  of  this  form,  one  has  /»(!->  )  =  qx\  where  *'  is  the  value 
at  maximum  i.  Using  this  relation  to  reduce  the  number  of  unknown  parameters 
to  one.  one  can  use  a  simpler  graph  to  determine  p  and  q  [26].  Flowever.  the  cal¬ 
culated  values  of  p  and  q  are  then  subject  to  the  accuracy  of  y  and.  more  funda¬ 
mentally.  the  possibility  that  p  and  q  may  change  from  one  range  of  y  to  another 
is  not  allowed  for.  As  mentioned  above,  those  types  of  fd  -  y)  in  which  q  -  Q 
represent  the  acceleratory  period  which,  in  general  should  be  followed  by  a  decay 
period  governed  by  a  different  form  of  fd  -  y). 


The  reaction  is  either  phase-boundary  controlled  or  unimolecular.  and  a  = 
=  K{\  ~  *)*.  as  seen  in  Table  1.  One  then  draws  the  graph  of  log  i  against 
log  (1  -  i)  to  find  r.  the  apparent  reaction  order. 

3.  w  i  0..*! 

The  reaction  is  diffusion  controlled  (see  Table  1).  One  has  to  test  separately 
whether  z(r)  is  parabolic  (the  diffusion  is  in  one  dimension),  or  —  i  Ind  -  y)  = 
=  A(l  -  xf  with  .f  =  0  (two  dimensions)  or  j  =  1  3  (three  dimensions). 

The  correlation  coefficient  in  the  least-squares  fit  serves,  by  measuring  the  li¬ 
nearity  of  the  y  graph,  as  a  quantitative  indication  of  the  confidence  to  be  attached 
to  the  identified  form  of /(I  —  xl.  It  may  be  that  p  and  q.  r.  or  s  change  once  or 
twice  as  the  reaction  proceeds  from  beginning  to  completion,  but  the  i  graphs 
will  show  it  by  displaying  several  linear  segments.  If  however,  a  part  of  the  graph 
say  from  (x,.  t,)  to  (x,.  t.,)  is  non-linear,  then  f(»)  has  changed  to  a  form  in  another 
group.  The  first  step  should  then  be  repeated  for  that  part;  log  [  — ln(l  -  [5:  -i-  a,]' 
[I  -  a,))]  is  plotted  1-5.  log(r  -  /,)  for  a  between  a,  -1-  O.I5(a...  -  a,)  and  a,  +  0.5 
(a.,,  -  a,),  followed  by  one  of  the  above  three  alternative  procedures.  On  the  other 
hand,  if  a  good  fit  is  found  with  values  of  p  and  q.  r.  or  .r  that  are  not  in  Table  1. 
the  experimenter  should  assess  whether  theoretical  justification  can  be  provided. 
In  this  way  new-  rate  laws  may  be  identified. 

Confirmation  is  carried  out  by  plotting  the  selected  functional  form  or  forms 
on  top  of  the  experimental  curve.  A  slight  misfit  in  the  very  early  part  (-  1 
minute,  depending  on  the  sample  .•>ize  and  the  environment)  may  be  attributed 
to  thermal  lag-time  and  ignored.  The  i  graphs  give  A",  and  from  a  set  of  A  values 
at  different  temperatures  A' ,  and  £  can  be  determined.  Note  that  the  determination 
of  these  macro-kinetic  constants  depends  on  the  form  of  F(  I  -  x)  chosen,  as  it 
should  be. 
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It  may  be  added  that  we  see  it  an  immediate  possibility  to  have  full  automation 
in  the  acquisition  and  processing  of  data  in  thermal  analysis  experiments.  The 
hardware  can  be  under  the  control  of  microprocessors  or  dedicated  minicomputers, 
and  their  output  would  go  into  a  computer  or  the  same  minicomputer.  A  com¬ 
puter  program  can  then  reduce  the  data  to  x(t)  or  i(r)  curves,  and  further  analyse 
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Fig.  1.  Decomposition  of  a  single  crystal  of  AgN,  at  551  K.;  1.  Thermogravimetnc  data. 
2.  a-timc  curve.  3.  1  -  (1  —  »}'  •  w.  time. 
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Fig.  2.  Log-In  plot 


the  curves  to  idcntify/^1  -  a)  and  so  calculate  and  £.  according  to  the  method 
proposed  here.  Nevertheless,  the  physical  interpretation  of  these  results  by  the 
experimenter  remains  the  crucial  step. 

The  method  has  been  applied  to  investigate  the  kinetics  of  slow  thermal  de¬ 
composition  in  silver  azide  single  crystals  using  thermogravimetric  data.  Curve  1 
in  Fig.  1  is  a  typical  experimental  curve  of  weight  loss  against  time  r,  and  Curve  2 
is  the  corresponding  reduced-time  plot  of  x  cs.  t  t.  In  Fig.  2  we  plot  logt— ln(l  - 
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—  2)]  w.  logU/f)  for  0.15  <  a  <  0.5.  obtaining  essentially  a  linear  graph  of  slope 
~1.1.  The  decomposition  of  AgN:,  therefore  appears  to  be  phase-boundary 
controlled  or  unimolecular.  Accordingly,  we  draw  in  Kig.  'S  the  graph  of  log  i  vs. 
logll  —  x).  It  show's  that  between  I  —  a  =  I  to  0.1.  i  =  K{\  —  a)-  with  the 
correlation  coefficient  among  the  data  points  better  than  -1-0.9.  Indeed,  the  plot 
of  1  —  (1  —  a)*  vs.  r  (Curve  3.  Fig.  1)  is  a  good  straight  line,  with  a  correlation 
coefficient  of  -(-0,99  for  0. 1  <  :)(  <  0.9.  The  decomposition  is  thus  of  the  contract¬ 
ing  cylinder  type.  We  have  also  obtained  K,  and  E.  and  found  that  their  inter¬ 
pretation  in  terms  of  Eqs  (3)  and  (4)  leads  to  a  plausible  physical  picture.  Further 
details  are  given  in  a  separate  paper  devoted  to  the  mechanism  of  the  thermal  de¬ 
composition  of  AgN:,  (^O)- 

In  a  second  paper  [29],  we  discuss  the  analysis  of  dynamic  kinetic  data  which 
are  easily  obtained  from  modern  thermoanalytical  instruments. 


Conclusion 

We  have  described  a  systematic  method  for  determining  the  kinetics  of  solid- 
state  reactions  from  isothermal  data.  It  may  be  stressed  once  more  that  whenever 
possible  a  judgement  should  then  be  made  on  the  consistency  of  the  AT®  gnd  E 
values  with  the  implication  of  /(I  —  a)  regarding  the  likely  mechanism  of  the 
reaction. 
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The  kinetics  of  heterogeneous  reactions,  involving  one  reactant  in  the  solid  phase, 
usMiUy  follow  the  law  a  =  ATo,  exp(— f/ikTT/Il  —  s),  where  a  is  the  degree  of  conversion 
of  the  solid,  and  and  E  are  the  kinetic  constants.  A  critical  examination  is  given  of 
(he  various  methods  which  are  currently  used  to  analyse  dynamic  experimental  data. 
The  limitations  of  these  methods  and  their  insensitivity  to  the  form  of  /(I  —  a)  are 
pointed  out.  An  alternative  approach  free  from  these  limitations  is  suggested.  In  this, 
/(I  —  a)  is  determined  from  isothermal  experiments,  and  then  the  dynamic  data  are 
accurately  analyzed  to  obtain  the  values  of  the  kinetic  constants.  A  case  study  is  given 
to  elucidate  the  applicability  of  the  approach. 

There  are  many  reactions  of  interest  in  which  one  of  the  reactants  is  in  the  solid 
phase.  These  reactions  can  be  classified  variously  as  decomposition,  dehydration, 
calcination,  dehydroxylation,  reduction,  polymeric  inversion  and  degradation, 
oxidation  etc.,  and  they  occur  in  a  wide  range  of  substances  including  ceramics, 
explosives  and  biological  materials.  The  dynamic  method  of  studying  their  kinetics 
involves  measuring  the  reaction  rates  under  conditions  of  a  continuous  temperature 
change. 

Originally  proposed  by  Skramovsky  [1  ],  the  dynamic  method  is  becoming  in¬ 
creasingly  popular,  especially  with  the  development  of  differential  thermoanalytical 
techniques  like  DSC  and  DTA.  Since  the  initial  temperature  can  be  chosen  so  that 
the  reaction  rate  is  relatively  insignificant  to  begin  with,  it  does  not  suffer  from 
zero-time  inaccuracy -a  problem  which  exists  in  isothermal  experiments  where  the 
temperature  is  raised  rapidly  and  then  held  constant  at  a  particular  value.  A  further 
advantage  is  that,  provided  the  dynamic  data  have  been  unambiguously  and  cor¬ 
rectly  analyzed,  any  changes  in  the  kinetic  constants  will  not  be  over-looked  even 
within  small  temperature  intervals.  In  contrast,  the  isothermal  method  only  pro¬ 
vides  values  averaged  over  discrete  points  in  temperature.  Also,  when  the  method  of 
analysis  used  is  such  that  the  kinetic  constants  are  calculated  from  each  dynamic 
curve  then  very  few  samples  are  required;  only  a  milligram  or  so  of  the  material  is 


needed  for  its  thermal  characterization.  If  many  runs  are  indeed  carried  out.  differ¬ 
ences  between  individual  samples  can  be  determined.  The  last  two  advantages  are 
particularly  useful  in  single  crystal  work. 

On  the  other  hand,  intrinsic  differences  should  be  carefully  distinguished  from 
the  effects  of  experimental  conditions.  In  the  first  case  literature  data  have  shown 

that  experimental  parameters  such  as  sample  mass  [S3]  and  shape,  particle  size  in 
the  case  of  powder  samples,  and  ambient  atmosphere  (or  vacuumican  affect  signif¬ 
icantly  the  calculated  values  of  kinetic  constants  (see  [2]).  If  this  happens,  then 
whenever  p<sssible  the  empirical  results  should  be  extrapolated  to  refer  to  a  stan- 
ilanl' set  ol'cxjscrimenial  conditions.  Secondly,  the  heating  rate  very  often  affects 
the  shape  of  the  dynamic  curve  obtained,  but  discussion  will  be  deferred  to  the  next 
section.  Lastly,  the  very  fact  that  temperature  is  now  a  variable,  in  addition  to  time, 
complicates  the  analysis  of  data.  If  due  care  is  not  taken,  either  inaccurate  or  totally 
misleading  values  are  obtained.  In  fact,  a  survey  of  the  literature  reveals  several 
instances  of  high-quality  experimental  data  being  mis-interpreted  by  methods 
beyond  their  ranges  of  validity,  in  this  paper,  we  first  describe  the  various  methods 
and  point  out  their  limitations.  We  then  defend  the  approach  in  which  use  is  made 
of  both  isothermal  and  dynamic  experiments.  The  analysis  of  isothermal  data  yields 
/(I  -x)  (as  defined  below)  unambiguously  and.  knowing/fl  -x).  one  can  calculate 
individual  values  of  the  kinetic  constants  from  each  set  of  dynamic  data.  We  may 
mention  that,  historically,  isothermal  experiments  were  the  only  ones  employed 
in  the  pioneering  age  in  the  twenties  and  thirties,  when  solid-state  reactions  began 
to  be  studied  from  the  modern  point  of  view,  as  distinct  from  that  of  Langmuir. 
Nernsi  and  Tammann. 

Kinetic  equation 

As  discussed  elsewhere  [3]  the  kinetics  of  a  reaction  proceeding  isothermally  can 
usually  be  described  by  the  empirical  relation ; 

(a)i«,ti«rm.i  =  Af,/(l-a)exp(-E/<:r).  (I')- 

Here  x  is  the  fraction  of  the  solid  reacted.  <c  Boltzmann's  constant.  T  the  tempera¬ 
ture.  and  /(I  -x)  and  the  constants  AT,  and  £  are  characteristic  to  the  reaction. 
The  function  /( 1  -  x)  may  change  in  different  ranges  of  x  but  is.  for  a  given  x.  inde¬ 
pendent  of  r.  at  least  within  a  range  of  T.  and  £  should  be  the  same  for  the 
same /(I  -x).  If  the  rate-controlling  step  of  the  reaction  occurs  on  the  reactant -free 
surface  or  on  the  reactant-product  (solid)  interface,  then  will  contain  the  sur- 
face-to-volume  ratio.  In  other  words,  the  reacting  system  should  really  be  normal¬ 
ized  per  unit  area  rather  than  per  unit  size,  and  be  given  in  units  such  as  mole¬ 
cules  s'’ m"'. 

Some  authors  have  questioned  the  general  validity  of  (1)  on  various  grounds 
[4-7].  However,  in  the  literature  (I)  is  almost  always  successfully  fitted  to  experi¬ 
mental  data.  Indeed,  this  empirical  relation  can  be  given  mechanismic  jutification 
(see  [3]). 

In  dynamic  experiments,  also,  it  is  commonly  agreed  that  (1 )  may  be  adapted  to 
describe  the  reaction  rate; 

X  =  Alt  /(1 -x)exp  [-f/A:7'(r)]  (2) 

£(()  is  controllable  by  the  experimenter.  Some  temperature  programs  offer  the 
mathematical  advantage  that  exp  {-£MT]/t  can  be  integrated  analytically  (see 
later  discussion  on  integral  methods  of  data  analysis).  Examples  are  the  hyperbolic 
program  where  LT  =  /(  -  fft.  i.c.  t  =  BT'  (8,  9];  a  parabt>lic  program  in  which 
AT~  +  BT  -  C  =  t.  i.e.  \iT  -  2AT  +  B.  with  B  =  AEiK,  by  iteration  [10]:  and 
an  exponential  program  so  that  t  =  exp  (-BIT)  with  B  &  E,k  [10].  For  the  sake 
of  experimental  convenience,  however,  the  arrangement  is  usually  that  f  »  ^ 
(A.  C  and  represent  constants  in  a  particular  run  of  experiment). 
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However,  it  has  been  taken  by  some  authors  who  object  to  (2).  that 
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where  ..  s  ( “)i,oihCTm»i-  The  argument  is  then  that  (2)  is  seen  to  be  inade- 

quate  even  in  the  case  of  d<f>  =  0.  since  the  second  term  on  the  R.H.S.  of  (3)  is  non¬ 
zero  but  left  out  in  (2j;  see  e.g.  (II.  I2j.  It  has  further  been  proposed  (13).  by  a 
derisaiion  starting  from  (3).  that  (2)  is  correct  only  if  it  includes  the  extra  factor 

(1  +  (I  -  TJT)ElkT] 

in  which  T„  =  r(fl).  But  we  hold  that  (3)  is  unsound.  Given  i.  T  and  0.  i  is  not 
uniquely  determined  and  therefore  not  a  function  of  these  system  variables:  0dx  > 
>  0.  as  a  cannot  decrease  even  for  negative  dT  and  d0!  Nevertheless,  the  inexact 
differential  dx  can  be  integrated,  if  the  dynamic  process  can  be  treated  as  the  limit¬ 
ing  case  of  a  series  of  time  intervals,  during  which  the  reaction  proceeds  isother- 
mally  according  to  (1)  but  at  the  end  of  each  of  which  T  is  altered,  in  a  time  so  short 
that  during  it  the  sample  is  unchanged.  Along  this  path  P  the  result  is  easily  ob¬ 
tained  [14]: 

,  T 

j.xp,-£,iridr  ,4, 

OPT, 

where  K  =  Kr  ex}p(-ElkT).  We  must  emphasize  that  (4)  is  not  logically  self- 
evident.  as  is  sometimes  implied  (14)  or  argued  by  mathematical  operations  based 
on  the  presumption  that  a  »  »(r.  t)  (15).  Rather,  it  comes  from  the  assumption 
that  the  reaction  under  study  involves  no  slow  processes,  so  that  i  depends  only 
on  the  present  values  of  a  and  r(ii  is  a  function  of  state),  but  not  on  the  history  of 
the  reacting  system  (c.f.  [16]).  Only  by  this  assumption  (absence  of  memory  effects) 
can  the  dynamic  process  be  treated  as  P.  (Experimentally,  a  temperature  program 
with  temperature  jumps,  which  approximates  P,  has  been  realized  on  a  thermo- 
balance  interactively  controlled  by  computer  [17].)  Equations  (2)  and  (4)  are  of 
course  equivalent.  Their  validity  has  also  been  shown  by  'rationar  thermodynamic 
arguments,  in  which  the  functional  relation  *  =  ,9  (t.  fit)  is  regarded  as  the  ‘con¬ 
stitutive  equation'  characterizing  the  reaction  system  [18].  In  a  new  direction. 

the  possibility  has  been  suggested  [64]  that  solid-state  reactions  may  be  studied 
by  far  from  equilibrium  thermodinamics. 

On  the  other  hand,  experience  shows  that  apparently  (4)  is  not  always  followed 
exactly.  Consider  a  reaction  being  investigated  by  a  series  of  experiments  conducted 
at  different  heating  rates  0  but  with  the  same  initial  temperature  T„.  In  (4)  we  see 
that  the  R.H.S..  for  a  given  upper  temperature  limit  T.  is  directly  proportional  to 
1/0.  Plots  of  the  L.H.S.  rs.  T  should  therefore  all  have  the  same  shape.  It  may  hap¬ 
pen.  however,  that  increasing  departure  from  isomorphism  is  seen  when  experimen¬ 
tal  data  obtained  at  higher  0  are  so  analyzed.  The  most  probable  explanation  is  that 
the  temperature  change  is  too  fast,  causing  the  temperature  distribution  in  the 
sample  to  become  significantly  non-uniform.  In  fact,  thermal  equilibrium  is  an 
underlying  assumption  when  (4)  is  derived  above:  without  it  i  will  depend  on  the 
thermal  history  of  the  reacting  system. 

Other  factors  may  also  be  at  work.  The  reaction  rate  may  be  sensitive  to  the 
structure  of  the  reactant,  and  a  higher  0  can  enhance  the  defect  density  or  change 
the  activation  energy  of  reaction  at  a  defect  site  (I9J.  In  branched-chain  reactions, 
the  speed  of  the  progressive  accumulation  of  active  centres  may  vary  with  the  rate 
of  change  in  temperature  [20].  If  the  reaction  is  a  surface  process,  the  distribution 
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ut  reaction  'centres'  among  comers,  edges  or  faces  of  the  sample  may  change  with 
<t>  [21].  It  may  also  be  that  the  chemical  system  under  study  has  multiple  reactions 
proceeding  concurrently  in  it,  and  they  have  different  £  [22].  All  these  variations 
in  £  may  be  accompanied  by  changes  of  in  the  same  direction.  Because  of  this 
coupling,  a  linear  relation  between  £  and  log  AT,  is  sometimes  observed.  Called  the 
‘compensation  effect’,  this  phenomena  does  not  necessarily  mean,  as  was  suggested 
[4],  that  the  Arrhenius  expression  in  the  R.H.S.  of  (4)  is  invalid.  In  all  these  cases, 
by  varying  0  the  experintenter  can.  in  fact,  gain  additional  insights  into  the  mecha- 
ni.sm  of  the  reaction,  or  distinguish  between  the  competitive  reactions  in  the  react¬ 
ing  systems  (u  situation  usually,  though  not  always,  indicated  by  the  presence  of 
multiple  peaks  in  the  i  curves).  This  is  possible  if  the  method  of  data  analysis 
employed  is  such  that  and  £  are  determined  from  a  single  i  or  a  curve,  rather 
than  from  data  at  a  number  of  healing  rules.  The  method  we  suggest  will  be  of  this 
type. 

It  may  also  happen  that  in  (2)  AT.  x  T,  so  that  it  cannot  be  taken  outside  the 
integral  sign  in  (4).  Indeed,  modifications  have  been  suggested  of  some  methods  of 
data  analysis  (those  that  assume  a  reaction  order  for  the  reaction)  to  take  this  extra 
temperature  dependence  into  consideration  [23].  However,  even  wlien  theoretically 
required,  the  correction  may  for  practical  purposes  be  ignored,  unless  £  is  small  or 
temperatures  used  are  very  high  (d  In  KjdT  «  [£  -h  kTykJ*).  Likewise,  any  slight 
temperature  dependence  of  £  can  usually  be  neglected.  Furthermore,  irrespective 
of  this  or  the  above  complications  the  form  of/f  1  -  a)  in  (2)  and  (4/  is  not  affected. 
It  should  be  the  same  as  that  in  (I).  on  the  basis  that  the  dynamic  process  can  be 
treated  as  the  limiting  case  of  a  series  of  isothermal  intervals,  as  already  mentioned 
above. 


Data  analysis 

Experimentally,  dl  or  a  is  obteined  by  DSC.  OTA,  TO,  DTG,  quantitative  IR 
spectroscopy  or  X-ray  diffraction,  dilatometry,  or  measurements  of  chemi-lumines- 
cence.  ultrasonic  attenuation,  dielectric  constant,  viscoelasticity,  thermal  or  electri¬ 
cal  conductivities,  or  optical  reflectivity  when  changes  in  these  characteristics  can 
be  correlated  with  x.  The  oldest  technique  is  thermomanometry,  in  which  the  pres¬ 
sure  of  an  evolved  gas  is  measured  at  constant  volume,  but  its  use  has  so  far  been 
more  popular  in  isothermal  experiments.  Many  methods  of  analyzing  d  or  x  data 
have  been  proposed  to  calculate  the  kinetic  consunts  £  and  AT.  (for  a  critical 
review  of  the  earlier  work  see  [24])  and  sometimes  also /(I  -  a).  Often  they  were 
originally  formulated  with  reference  to  one  particular  instrumentation,  but  they 
may  be  made  generally  applicable  to  all  techniques  afler  quantities  measured  on 
DSC,  DTA.  TG  instruments  etc.  are  all  interpreted  in  terms  of  xaad  i.  On  the 
other  hand,  their  validity  does  depend  on  the  particular  reaction  whose  data  ate 
being  analyzed.  Their  limitations  in  this  respect  form  the  subject  of  our  discussion 
below.  They  will  be  examined  in  three  groups:  peak-temperature,  integral  and 
derivative  methods,  in  this  order.  Sophisticated  instrumentation  systems  are  com¬ 
ing  into  use,  that  incorporate  computers  to  establish  baselines  or  other  null  settings, 
to  carry  out  automatic  data  acquisition,  and  to  let  the  experimenter  interactively 
analyse  the  data  (e.g.  [23]).  Such  advances  do  not,  however,  remove  the  danger  of 
uncritical  choices  of  the  method  of  data  reduction. 
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Critical  examlaatioB  of  cnrreu  methods 

Peak-temperature  method 

Kissinger  [26]  considers  reactions  of  the  type/Kl  -a)  *  (1  -a)*.  Differentiating 
(2)  with  respect  to  t,  and  setting  the  resulting  expression  to  zero,  he  obtains 

iJElkUITl)  -  txp(-ElkT„)n(l  -a,)*  '**  (5) 

in  which  m.  signifies  ‘peak*  quantities,  at  the  point  of  maximum  d  where  a  =  0. 
He  next  that  /»(!  —  a„)"  '  a:  I ;  therefore 

tftTjoLexpi-ElkTJ  (6) 

regardless  of  /?,  which  itself  may  be  calculated  from  the  shape  of  the  a  (t)  curve. 
E  anil  k.  on  the  other  hand,  are  obtained  by  performing  a  series  of  experiments  at 
different  0.  An  aspect,  which  we  regard  as  an  inefficiency,  of  Kissinger's  method  is 
that  only  one  point  on  the  curve  is  used  although,  in  the  case  where  multiple  peaks 
occur  signifying  that  different /( I  —  a)  and  E  govern  different  sections  of  the  curve, 
the  method  should  still  be  applicable  to  each  peak. 

There  is,  however,  an  important  limitation.  The  a  priori  condition  that/fl  -  x)  > 
s  (I  -  a)”  is  actually  valid  only  in  very  special  circumstances,  namely  when  the 
rate-limiting  step  of  the  reaction  is  the  inward  movement  at  a  constant  speed  of 
the  reactant-product  interface,  where  n  is  0.  V->  or  2/3  for  one-,  two-,  or  three- 
dimensional  movement,  respectively,  or  when  the  reaction  is  unimolecular  so  that 
n  =  [.  Even  among  these  special  cases,  the  other  approximation  that  Kissinger 
uses  is  still  conditional,  since  n(l  -a„)"' '  s  1  only  for  «  =  1.  When  n  is  V:  or  2/3. 
this  expression  varies  with  approximately  as  da„,  n(l -ti)/(l *  0.2 

da„.  where  is  the  variation  in  itself.  In  the  Appendix  we  show  that  a„ 
changes  with  0  in  the  general  case.  Hence,  when  an  apparent  reaction  order  exists 
and  is  '/i  or  2/3,  Kissinger's  method  can  lead  to  a  systematic  deviation  in  (6)  and  thus 
generate  a  significant  but  hidden  error  in  the  calculated  Eand  K, . 

If  no  apparent  reaction  order  exists,  then  it  definitely  should  not  be  used,  other¬ 
wise  an  approximately  linear  plot  from  (6)  results  in  totally  misleading  values  of 
the  kinetic  constants.  An  example  is  in  the  decomposition  of  benzenediazonium 
chloride;  it  derives  from  DTA  data  a  value  off  that  is  40%  lower  than  the  nearly 
identical  values,  obtained  by  applying  other  methods  of  analysis  to  the  data  from 
DTA  as  well  as  other  techniques  [27J.  Other  examples  are  in  the  study  of  lithium 
aluminium  hydride,  where  the  Kissinger  values  are  half  of  the  isothermal  result 
[28],  in  ROX  where  it  is  again  40%  lower  than  all  the  values  calculated  by  other 
methods  [29],  and  in  urea  nitrate,  where  it  is  30%  lower  [30]. 

Integral  methods 

The  L.H.S.  of  (4)  is  a  function  of  a  only  and  will  be  denoted  by  F(a):  the  R.H.S. 

r 

can  for  practical  purposes  be  equated  with  [  X/0  dT,  since  in  experiments  T„  will 

b 

be  such  that  reaction  velocity  is  negligible  below  it.  i.e.  T„  <<  Ejk.  In  view  of  these 
considerations,  many  authors  have  proposed  different  methods  of  analysing  aUT) 
data.  r 

The  temperature  integral  )  exp  {—ElkTfiT  has  no  analytical  solution.  (In  the 
b 

unusual  case  of  a  hyperbolic,  parabolic  or  exponential  temperature  program, 
on  the  other  hand,  exp  (-ElkT)/t  is  integrable.)The  numerical  values  of  the  inte¬ 
gral  have  been  compiled  but,  being  a  function  of  both  E  and  T,  arc  not  directly 
useful  unless  an  iterative  solution  of  (4)  by  trial-and-error  is  resorted  to.  Such  an 

approach  has  been  advocated  by  Zsakd  [31  ]  who  considers  in  particular  the  cases 

% 

of  .A  I  -  a)  =  (I  —  at)*  with  w  —  0,  !/3,  1/2,  2/3,  I  or  2,  when  log  j  — ■■ 

./  /I  I  —  «) 

u 
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has  simple  analytical  expreawons,  and  by  Satava  and  Skv&ra,  [32]  who  generaHze 
the  method  slightly  by  tabulating  the  values  of  p(a).  0  <  a  <  I,  for  some  other 
forms  of/]  I  -  a). 

For  more  efficient  approaches,  approximations  to  the  integral  are  necessary. 
Thus,  talcing  the  first  two  terms  in  an  asymptotic  (u  s  E/kT  -»  oo)  series  of 

k 

r  dr 

I  exp  ( -  t)  Coats  and  Redfern  obtain  the  linearized  relation  [33] : 

In  (F]a)7"  ^)  =  .4  -  £/*(r ' ' )  (7) 

wheic  A  =  ln(A.’ ,  \  —  24  77£') is  "sensibly  constant"  if  the  range  of  temperature 

AT  is  small.  They  further  assume  that /Hi  -a)  »  ( I  -a)*,  and  so  f(x)  can  be  cal¬ 
culated  at  each  (a,T).  Plotting  (7)  for  several  values  of  a  thus  gives  E  and  AC,, 
making  use  of  only  one  set  of  data  ct>rrcsponding  to  a  single 
Several  cautionary  notes  should  again  be  made  here.  The  assumption  for/]!  -  a) 
has  already  been  discussed.  Similar  to  the  case  of  Kissinger's  method,  resulu  ob¬ 
tained  may  be  wrong  and  misleading  if  this  functional  form  is  not  independently 
determined  beforehand.  Thus,  in  a  study  on  the  dehydroxylation  of  kaolinite  [34], 
straight  lines  over  different  ranges  of  (T-*)  are  given  by  (7)  for  a  whole  series  <rf 
values  of  n.  namely,  0,  O.S.  0.667.  a  and  2.  In  particular,  plots  using  n  *  I  and 
n  =  2  are  almost  equally ‘good*. 

Secondly,  the  accuracy  of  the  asymptotic  approximation  is  rather  low.  By  com¬ 
paring  its  values  with  ttbulated  values  of  the  integral  (31.  35.  36],  we  find  its  rela¬ 
tive  errors  to  be  4///  -  20 %  at  u  -  5,  5%  at  u  -  10.  and  1 .5  %  at  u  -  20.  Thus,  for 
example,  if  £  is  I  eV,  then  for  an  accuracy  of  98  %  the  highest  temperature  reached 
in  the  experimental  run  should  not  be  more  than  600  K.  a  very  low  figure  for  most 
materials  though  it  is  higher  for  larger  £.  Additionally,  expanding  A  into  a  power 
series  shows  that44/4a!2/r477£,so  that  at  say  2?^  inaccuracy  the  range  of  tempera¬ 
ture.  AT.  from  which  (a,  T)  points  arc  selected  should  be  less  than  lOOK  (for  £  « 

1  eV).  The  total  possible  deviations  in  the  calculated  £  and  ,  are.  to  first  approxi¬ 
mation.  the  sum  of  the  4/// and  44/4.  It  certainly  is  unsatisfactory  if  they  ate  hirgd 
and  yet  nowhere  mentioned  in  the  calculation. 

Other  approximations  to  the  temperature  integral  have  been  suggested  by  van 
Krevelen  er  at.  [37]  and  by  Horowitz  and  Metzger  (38).  who  made  use  of  certain 
asymptotic  expansions  in  the  vicinity  of  T*.  the  temperature  at  peak  reaction  rate. 
Both  have  been  shown  (39)  to  be  even  less  accurate  than  the  Coats  and  Redfern 
approach,  and  so  will  be  left  out  in  our  discussion. 

Amongst  the  integral  methods,  the  best  is  probably  the  one  due  to  Ozawa,  which 
requires  data  at  different  ^  but.  in  it/]l  —  z)  remains  completely  general.  The  ap¬ 
proximation  to  the  temperature  integral  is :  - 

r 

J  exp( -£/*Ddr a IC 

SO  that  from  (4) 

log  ♦,  +  0.457(£/*)/r,  -  log  +  0.4S7(£/*)r,  (9) 

where  T,  and  7|  arc  uken  at  an  arbitrary  but  identical  value  of  a  in  the  two  curves 
corresponding  to  heating  rales  and  Plotting  log  ^  ps.  1/7"  for  selected  values 
of  a  should  therefore  produce  straight  lines,  the  slopes  of  which  give  £  (40). 

Three  comments  are  appropriate  here.  By  comparing  (8)  with  ubuhled  numer^ 
ical  values,  we  see  that  it  is  7%  out  at  w  -  10 or  7"  -  I I70K,  and  3%  and  leu  only 
for  T  <  720K  (if  E  a  1  eV).  These  errors  should  be  examined  before  Ozawa’s 
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method  is  applied.  Secondly,  the  method  has  been  modified  [41]  to  read,  in  place 
of  (9). 

d«n(^/dr„'  -  0.457  £/)i  (10) 

in  which  J  In  d*  ^  In  (t>i  —  In  <^>3,  etc.,  and  m  denotes,  as  before,  peak  quantities. 
This  relation  may  be  compared  with  (6)  but  in  general  it  does  not  hold  since,  as 
shown  in  the  Appendix.  varies  with  Lastly,  like  Kissinger's  method.  £  cannot 
be  determined  from  data  at  a  single  tp.  and  in  some  cases  this  may  be  adisadvantage. 
as  discussed  before. 

Derivative  methods 

The  derivative  methods  offer  an  advantage  over  those  described  above  in  invok¬ 
ing  no  mathematical  approximations.  Unfortunately,  they  use  i  data  which,  with 
present  instrumentation,  tend  to  be  of  lower  quality  whether  they  are  obtained  by 
numerically  dilfercniiating  the  1  data  or  are  direct  experimental  read-outs. 

The  most  straightforward,  but  as  it  stands  relatively  inefficient,  of  the  derivative 
methods  is  to  write  (2)  as;  — 

ln(d//(l-a))  =  In  K,  -  EIkT  (II) 

and  to  substitute  different  of/(  I  -  x)  until  a  linear  plot  appears  [42].  Later,  we  shall 
argue,  however,  that  even  this  labor  omnia  vincit  approach  like  all  dynamic  meth¬ 
ods  in  general,  cannot  guarantee  correct  values  of  £  and  K,  (nor  an  unambig¬ 
uous  form  of  /i[l  -x)  in  this  specific  case),  although  the  labour  it  involves  may  be 
undertaken  by  the  computer. 

The  earliest  derivative  method  is  probably  that  of  Borchardt  and  Daniels,  origi¬ 
nally  formulated  for  homogeneous  reactions  in  the  liquid  phase  [43]  but  later  extend¬ 
ed  to  solid-state  reactions  [44]  for  which  it  is  now  frequently  used.  The  method 
puts/]  I  -  X)  »  ( I  -  x)”  into  (ID,  with  n  given  a  guessed  value,  and  if  a  linear  plot 
results  then  £  and  Ky  are  obtained  from  it.  Based  on  this  method.  Hauser  and 
Field  [45]  have  developed  a  computer  procedure,  in  which  plots  are  generated  for 
a  series  of  values  of  n  incremented  at  discrete  steps,  and  the  'best'  one  is  then  se¬ 
lected  to  yield  £.  Ky.  and  n.  An  attraction  of  this  method  is  that  n  can  be  readily 
selected  by  eye.  Alternatively,  since  in  this  case 

J  In  X  dln(l -X)  »  — (£  A[dT' ' d  In  (I -x)]) -I- n  (12) 

a  plot  of  the  L.H.S.  vs.  the  quantity  in  the  square  brackets  at  once  gives  £  from  the 
slope  and  n  as  the  r-intercept[47].  If  constant  dlnx.  dln(l  -x).or  dT' '  is  selected. 
Eq.  (12)  can  be  further  simplified  [61).  We  have  emphasized  previously  the  falli¬ 
bility  in  presuming  such  a  convenient  form  of /(I  — x):  Ozawa  [47]  has  commented 
on  the  possibility  that  this  procedure,  and  the  integral  method  of  Coats  and  Red- 
fern.  may  give  false  values  of  £  and  AT, .  In  addition,  since  (12)  involves  the  ratios 
of  differences,  the  quality  of  data  called  for  is  even  higher  than  that  demanded 
alone  by  the  use  of  x ;  experimental  data  so  plotted  more  often  than  not  show  very 
large  scatter.  The  Rogers  and  Morris  metht^  [48]  plots  dlnx  r.v.  T~ and  can  be 

seen  to  be  the  special  case  of  «  =  0  in  (12).  An  example  of  the  general  danger  that 
very  linear  plots  may  sometimes  appear  even  if  the  applied  method  is  not  valid  is 
given  by  Patel  and  Chaudhri.  The  Rogers  ar  d  Morris  method  was  used  to  analyze 
DSC  data  on  lead  azide,  and  a  straight  line  results  although  the  calculated  £  turns 
out  to  be  180“o  larger  than  the  Ozawa  value  [49].  Conversely,  the  coincidence  of 
values  calculated  by  various  methods  need  not  prove  that  these  methods  are  all 
applicable  to  the  case  in  hand.  A  counter-example  is  provided  by  a  DSC  study  on 
RDX  [29].  where  the  Rogers  and  Morris  value  agrees  well  with  other  values'but 
the  compDx  decomposition  is  beyond  doubt  far  from  the  n  ■  0  type. 

On  the  other  hand.  D4vid  and  Zelenyanszki  [50]  plot  In  j  (1  -  x)j/(l  -  x) 

against  (£*');  this  amounts  to  assuming  a  'reaction  order'  n  =  1.  It  serves  as  yet 
another  example  of  the  futility  of  linear  plots,  for  their  method  gives  such  plots  for 
the  decomposition  of  'a  wide  range  of  materials'  including  calcium  oxalate  and 
polyethylene  which,  most  likely,  are  not  of  first  or  any  other  'order'. 
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Some  of  the  integral  and  derivative  methods  described  in  the  foregoing  have  been 
compared  by  testing  their  accuracies  on  synthesised  DTA  data  (exact  as  well  as 
with  artificial  random  error)  for  one  E  value  and  temperature  range,  the  reaction 
considered  being  of  the  type  with  a  reaction  order  [51].  Among  the  methods  not 
included  there  is  that  due  to  Friedmann  [52].  It  probably  is  the  most  general  among 
the  derivative  methods.  Like  Ozawa's  procedure,  it  makes  no  assumption  about 
/( 1-2).  although  it  requires  a  data  which,  furthermore,  have  to  be  at  a  number  of 
Once  again,  from  (2)  with  da  dT"  =  a^: 

=  \n(K,f(l-ot))  -  ElkT.  (13) 

Since  A ,  I  —  a)  is  identical  for  the  same  value  of  a,  taking  and  the  correspond¬ 
ing  Ffrom  several  tj)  one  can  determine  E. 

Ii  is  our  conicmion  that  even  Friedmann's  method  has  one  basic  limitation  which, 
more  significantly,  is  shared  by  all  dynamic  methods  described  above.  The  point 
in  question  is  that  all  of  them  have  to  presume  the  constancy  off(  I  —  a)  as  the  tem¬ 
perature  IS  ehanged.  However,  since  mechanisms  of  solid-state  reactions  are  gener¬ 
ally  complicated,  there  is  no  general  justification  for  this  presumption,  though  it 
may  be  true  for  particular  reactions  within  specific  temperature  ranges.  An  illustra¬ 
tion  is  the  case  where  parallel  reaction  paths  exist,  each  with  values  of  and  E 
such  that  a  quantitative  change  of  T  will  lead  to  a  qualitative  change  in  the  dominat¬ 
ing  path.  Another  case  is  where  the  identification  of  the  rate-limiting  step  depends 
on  T.  Methods  have  been  proposed  which,  by  the  use  of  computers,  try  different 
forms  of  /(I  -x)  in  analyzing  the  dynamic  data  [53-55].  However,  the  search  is 
limited  to  functional  forms  which  are  already  known. 

More  importantly,  from  our  own  experience  with  azides  we  have  strong  doubts 
as  to  the  exactness  in  determining/(  1  —  a)  or  even  its  constancy  from  dynamic  data. 
Likewise,  in  a  study  on  the  dehydration  of  manganese  formate  [54]  for  instance, 
no  unique  form  of  /(I  -  a)  and  correspondingly  no  unique  values  off  are  identified 
even  over  appropriately  restricted  ranges  of  a,  the  criterion  used  being  minimum 
standard  deviation  in  the  Arrhenius  plot.  Further  examples  are  the  thermal  dehy- 
droxylations  of  kaolinite  [34]  and  of  magnesium  hydroxide  [56].  We  suggest  that, 
in  dynamic  experiments  since  data  are  collected  under  variable  temperature  condi¬ 
tions,  the  change  due  to  /(I  -  a)  is  inherently  masked  by  that  due  to>C(7').This  pitfall 
is  illustrated  in  the  Figures.  Figure  la  shows  the  graphs  of  a(T)  and  its  derivative 
which  are  generated  artifically  according  to  the  theoretical  equation  a*"  = 


10' /T,  K" 


Fig.  ta.  Artificial  data  t  »  K/.  and  corresponding  i''  e  dr/drdata,  plotted  against  T  which 

rises  linearly  with  f 

Fig.  lb.  Arrhenius  plots  of;  (I)  »'^/3(l  -  (2)  *^,'2(1  -  »)"*.  (3)  and  ',(4)  iV3»*/* 

for  the  data  shown  in  Fig.  la 
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exp  ( —EjkT)  throughout  from  a  =  0  to  a  =  1 .  The  values  chosen  for  and  E  are 
10*  and  leVrespectively.  Let  us  now  examine  how  the  data  generated  according  to 
this  relation,  which  is  of  the  type  a  =  Ki,  will  be  fitted  by  different  kinetic  equations, 
one  of  them  being  the  correct  one.  In  Fig.  lb  we  plot  against  10*/r  the  natural 
logarithms  of  the  following  expressions:  — 

(1)  aV3(l— a)*'*,  i.e.  assuming  1  —  (1— a)^'*=  Kt:  reaction  controlled  by 
three-dimensional  contraction  of  phase  boundary : 

(2)  t^/20  -  a)'  ",  i.e.  1  -  (1  -  a)'  ■  =  /ft  type: 

(3 1  a^:  the  original  assumption :  and 

(4)  a73a‘''‘,  i.e.  a  =  (Af/)'*:  reaction  controlled  by  e.g.  three-dimensional 

growth  of  existing  nuclei. 

It  is  seen  that  the  incorrect  /HI  —  a)  in  (I land  (2) still  give  virtually  linear  plots,  with 
slightly  different  slopes:  interestingly,  curve  (4)  is  so  misleading  as  to  show  two 
“linear"'  segments  with  a  seemingly  significant  transition  in  between.  Experimen¬ 
tally,  Guarini  etal.  [57]  have  noted  that  it  is  impossible  to  ascertain  from  their  DSC 
data  whether  the  monomerization  of  9-Me-IO-AcAD  has  an  apparent  reaction 
order  of  1 .  0.67.  or  0.5,  in  all  of  which  cases  £  has  about  the  same  derived  value. 

The  suggested  approach 

In  the  foregoing  sections,  we  have  discussed  the  limitations  regarding  the  applica¬ 
bility  of  various  methods  that  have  been  used  to  analyse  dynamic  data.  In  many 
published  works  we  find  that  often  a  number  of  apparently  different  methods  are 
used  to  analyse  the  same  set  of  data.  However,  we  think  that  in  many  cases  this 
procedure  is  of  no  real  significance,  when  some  of  the  methods  used  are  mathemat¬ 
ically  equivalent  and  therefore  lead  to  the  same  results,  or  when  some  are  invalid 
in  Che  given  situation  and  thus  lead  to  doubtful  values.  The  limitations  of  the  meth¬ 
ods  express  themselves  both  as  discrepancies  in  the  calculated  values  of  the  kinetic 
constants,  and  sometimes  as  fortuitous  agreements  when  some  of  the  methods  are 
certainly  inapplicable.  (An  extreme  example  of  the  second  situation  is  that,  for 
RDX.  the  Kissinger  value  [29]  of  £  is  near  to  that  obtained  [58]  by  plotting 
2  rs  T~'.  a  procedure  which  has  absolutely  no  theoretical  justification.)  Accord¬ 
ingly.  we  suggest  that  the  interpretation  of  dynamic  data  should  as  far  as  possible 
be  based  on  results  from  isothermal  experiments.  A  similar  approach  has  been 
used  for  studying  the  dehydroxylation  of  kaolinite  by  Achar.  Brindley  and  Sharp 
[63]. 

One  can  unambiguously  determine /( 1  -x)  over  the  whole  range  of  x  and  over 
the  relevant  temperature  range,  from  the  independent  analysis  of  individual  iso¬ 
thermal  curves.  A  systematic  method  of  efficiently  implementing  this  identification 
has  been  proposed  by  us  [3].  It  may  also  be  noted  that  thermoanalytical  equipments 
are  equally  applicable  in  isothermal  experiments  (see  e.g.  [62])  though  they  are 
more  often  used  in  the  dynamic  mode.  The  identified  fonn(r)  of/Kl  -x)  can  then 
be  substituted  into  either  (2)  or  (4).  In  this  way.  form  the  dynamic  x  or  d  data  one 
can  then  determine  accurately  the  non-average  and  single-sample  values  of  £  and 
Kr :  advantages  which  have  been  mentioned  in  the  introduction  to  this  paper. 
Moreover,  the  values  will  correspond  individually  to  different  heating  rates. 

We  applied  this  approach  to  the  spinel  formation  ZnO  -t-  CrjO,  -*  ZnCr204. 
A  DTA  curve  (experimental  atmosphere:  Nj  at  300  mm  mercury)  was  published 
in  Ishii  etal.  [59],  who  have  also  monitored  x(0  by  chemical  analysis  when  the 
reaction  proceeded  isothermally  in  nitrogen  flowing  at  50  ml/min,  and  showed 
that  the  isothermal  data  fit  [I  -  (1  -x)''*]'^  «  Kt.  We  have  measured  K  from  the 
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I&olhermal  data 

r. 

deg.  C 

K 

relative  unit. 

In  k 

800 

1 

0 

900 

4 

1.4 

1000 

23 

3.1 

Tabic  2 

Dynamic  data 

/. 

lill.X,  II 

deg.  t' 

rctaiivc  units 

<  II  al 

700 

0.1  2 

-2.6 

820 

O.IS  4 

-1.4 

900 

0.2  7 

-0.54 

1000 

0.35  12 

+0.76 

experimental  data  points  at  r  =  20  min  in  the  published  isothermal  plots.  From 
these  values  of  K,  given  in  Table  1.  we  calculate  a  value  of  1.5  eV  for  £. 

In  Table  2  the  values  i  were  measured  from  the  published  DTA  curve  whose 
heating  rate  was  unspecified,  and  the  <x  values  were  read  off  from  the  (a,  T)  graph 
which  Ishii  el  al.  have  drawn  presumably  by  integration.  Now,  from  their  analysis 
of  the  isothermal  data  the  governing  kinetic  equation  is,  in  differential  form,  i  - 
=  A/((l  -  a)"*'*  -  (1  -  a)'*'®),  at  least  within  the  ranges  800- 1000°  and  a  from 
0  to  i  0.6  corresponding  to  £/  =  0  to  &  0.07.  The  Arrhenius  plot  of  i  (1  -  a)"  *'® 
((I -a)"'  ®  -  II,  for  the  four  data  points  shown  in  Table  2,  is  indeed  a  good 
straight  line.  From  the  plot  we  obtain  £  ==  1.3  eV.  In  view  of  the  probable  experi¬ 
mental  errors  and  inaccuracies  in  obtaining  data  from  the  published  graphs,  we  con¬ 
sider  satisfactory  the  reasonable  agreement  between  this  value  and  the  one  calculat¬ 
ed  from  the  isothermal  data. 


Coachision 

Most  of  the  commonly  used  methods  of  analyzing  dynamic  data  have  been  shown 
to  be  applicable  only  under  particular  conditions.  It  has  been  pointed  out  that  to 
use  these  methods  without  considerations  of  the  range  of  their  validity  can  give 
misleading  values  of  the  kinetic  parameters.  An  approach  has  been  advocated  in 
which  use  is  made  of  both  the  dynamic  and  isothermal  data:  the  functional  form 
is  determined  from  the  isothermal  experiments.  This  form  in  conjunction 
with  the  dynamic  data,  gives  the  values  of  the  kinetic  constants. 
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Appendix 

We  wish  to  predict  how  the  value  of  a  at  peak  reaction  rate  varies  with  the  heat¬ 
ing  rate  0.  At  *  =  a„:  from  (4) 


dot 

J  /(T^ 


0  T. 

and  from  the  fact  that  3  =  0 


J  dTexp{-ElkT)s  F(a„.  <t>.  TJ  =  0 


fi\  ,  »  ^  ^  t%p(ElkT„) 

f  (I  *>«)  ^  y.;  —  G(ct„.  (I>,  T„)  —  0 . 

Solving  the  simultaneous  equations  dF  =  0  and  dC  =  0.  we  find 


(I4> 

(15) 


d3„,  /  I'F  i)G  rF 

i'G  j  (  cF  <x;  ,'iF  .t;  j 

d</>  "  <'T,„ 

''(il  l  f‘T„  i'x„  i)T„  j  ■ 

(16) 

Defining  dimensionless  quantities 

U  =  ElkT„.  n  =  <l>IK,T„.  and 

*m 

lixj  s  I 

d3//(l  —  3).  we  have  the  following; 

b 

cF  _  1 

cG 

(‘x„  /([I-3J  ’ 

V- - /'(1-aJ. 

^F 


<iTcxp(-U)  = 

<4>  <p 


cF 


To 

I 


/(3,„)  cG  HU 


cG  nV(l  -t^) 


cT„  ntxp(U)T„  '  ST„ 

Hence  — ^  -  I 

d0  <t>  hlj(i  +  V)m- »«)  ~/*(l  -»„)//» 


exp(t/).  (17> 
(18) 


where  h  =  n  exp  {U).  (ncidentally,  dTJd<t>  can  tc  derived  in  a  similar  way. 

The  only  case  we  find  reported  in  the  literature,  in  which  3„  is  apparently  inde¬ 
pendent  of  0.  is  the  primary  recrystallization  of  pre-cempressed  copper  [60],  where 
x„  s  0.5.  In  all  other  cases,  experiments  give  changing  r„.  We  have  made  a  rough 
check  on  ( 1 8)  by  taking  the  case  of  the  decomposition  of  the  explosive  RDX  [29], 
for  which  the  Rogeraand  Morris  method  eives  E  =  2.10  eV  and  K.,  »  10**’*j“*. 
The  reaction  is  complex,  but  these  representative  values  are  chosen  because  they 
correspond  to  an  assumed  kinetic  equation  in  which/l[l  -3)  =  I.  We  thus  have 
very  simply /(3„)  =  3„and/*(l-«„)  =  0.  For  ^  =  0.167  K  s■^  3„  is  given  as  0.62 
and  T„  as  5I2K;  our  calculation  shows  d!xjd<l>  21  -0.1  s  K~'.  a  value  which 
compares  well  with  the  experimental  indication  that  =  (O-fiO  -  0.62)/ 

; (0.333  -  0.167)sK''. 
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